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Concerning the Division of Electron and Ion Optics 
of the American Physical Society 


By J. R. PIERCE 


Secretary-Treasurer, Division of Electron and Ion Optics 


RIOR to the inaugural meeting of the new 
Division of Electron and lon Optics, held 
at Pittsburgh on April 29 of this year, the 
members of the Division were invited to express 
their particular interests in the field and to seek 
and to contribute to clarification of the aims of 
the Division. | believe that the questions raised 
may be of interest to non-members as well as 
present members of the Division. Several general 
questions were raised and these perhaps deserve 
first attention. 

Generally, the aims and purposes of the 
Division are stated in the by-laws to be “the 
advancement of knowledge of electron and ion 
optics. To achieve this, the Division shall arrange 
for the presentation, at meetings of the Society 
or meetings of its own, of papers and discussions 
lying within its scope.’’ Such activities may well 
be supplemented by encouragement of publica- 
tion and other dissemination of information 
through correspondence and other activity on 
the part of the membership and officers. The 
Division plans no special publications. When 
meetings are held concurrently with general 
meetings of the American Physical Society, 
abstracts will appear in the Bulletin and be 
reprinted in The Physical Review. 

The scope of the Division is such as to include 
both theoretical and practical aspects of pro- 


duction and utilization of ion and molecular 
beams. It should certainly include consideration 
of both low frequency and high frequency aspects 
of electron flow, whether that flow occurs in a 
wholly new device, a mass spectrograph, or a 
vacuum tube. Many members of the Division 
will be interested in new developments and 
techniques in the use of the electron microscope. 
However, the particular field which no other 
organization seems to include so directly is the 
internal physics of electronic and ionic devices. 

One member wants to know whether the aims 
and purposes of the Division are beneficial to 
the parent society. I certainly hope so. Another 
is concerned with the social import of this branch 
of physics, of which I can say little. One is 
interested in management of research and de- 
velopment, a very concrete and vitally important 
problem of no little difficulty. 

Turning to specific technical matters, by far 
the greatest number of replies indicated a pri- 
mary interest in what might be termed general 
electron and ion optics. Interest was expressed 
in the primarily mathematical aspects of electron 
optics and in mathematics from the point of 
view of design, in electric and magnetic lenses 
and focusing means, in relativistic electron optics, 
especially in connection with magnetic analyzers 
for beta-particle spectra, in electron and ion 
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sources, in beams with low energy spread, and 
in molecular beams. 

More specifically, the electron microscope and 
electron diffraction were very frequently men- 
tioned. Other members expressed particular in- 
terest in the mass spectrograph, the cyclotron, 
x-ray tubes, and the electron optics of radio 
tubes. 

One member showed a particular interest in 
oxide-coated cathodes. | do not know whether 
this puzzling subject, together with the general 
study of thermionic emission, is naturally within 
the scope of the Division; I believe that it well 
may be. It is certainly a field in which much 
remains to be done. 

As examples of how the activities of the 
Division may be of help to the members, | 
mention that on the program of the Pittsburgh 
meeting we had invited papers covering several 
of the subjects in which interest was expressed, 
including the historical background of the topic 
of most interest, general electron optics, and 
papers on the perhaps somewhat neglected optics 
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of vacuum tubes as well as papers on the very 
popular subject of electron microscopy. 

As another example of such service as might 
be rendered through the Division, | refer to an 
excellent bibliography of the electron microscope, 
which was recently prepared in another connec- 
tion by Claire Marton, the wife of the Vice 
Chairman of the Division, and Samuel Sass. | 
think it should be a duty and activity of the 
Division to encourage the preparation of such 
bibliographies, of summaries of fields such as 
have appeared in Reviews of Modern Physics and 
in the Journal of Applied Physics, and also of 
simple tutorial papers which may serve as 
introductions to the various aspects of electron 
and ion optics which I have mentioned earlier. 

Finally, one member has suggested that after 
the war a symposium might be held to consider 
new techniques and data that have not been 
released during the war. | certainly hope that 
such a symposium can be held soon. In the 
meantime, it may be possible to hold at a future 
meeting a symposium on some less secret but 
nonetheless vital matter. 
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New Books 








Mathematics of Modern Engineering, Volume II 


By Ernest G. KELLER. Pp. 309+xii, Figs. 87, 
154X24 cm. John Wiley & Sons, Inc., New York, 
1942. Price $4.00. 


This book, second in a series of three, considers discrete 
physical problems which lead to ordinary differential 
equations. It is stated that the material and methods 
included have evolved out of research and the advanced 
course in engineering at the General Electric Company. 

The principal purpose of the book is to provide by 
examples an introduction to the methods of mathematical 
engineering, this field being described as a new phase of 
engineering, and defined as those portions of engineering 
which can be formulated mathematically. As in mathe- 
matical physics, the fundamental method of mathematical 
engineering consists of: (1) the reduction of the physical 
phenomena involved to a mathematical system amenable 
to treatment; and (2) the solution of the system. Step 
one, usually the most difficult, requires originality, insight, 
and a broad knowledge of the physical and mathematical 
sciences. Step two, while often difficult, can usually be 
carried through by known mathematical technique. 

The value of this book lies in its collection of problems 
illustrating the application of mathematics to certain types 
of engineering problems. The examples relating to loco- 
motives reflect the author’s experience and are particularly 
interesting. 

The general theory is treated rather unsystematically. 
Considerable improvement would have resulted if the 
fundamental definitions and principles had been introduced 
accurately and early in the exposition, rather than (as in 
a number of instances) in the style of an afterthought. 
Thus: the various definitions of non-holonomic systems 
(pp. 20-29); the use of Koenig’s theorem before its state- 
ment (never adequately given, pp. 20-1, 31, 44); the 
treatment of vector time derivatives given is often inade- 
quate (even for approximations, pp. 45-7). The inclusion 
of Kron’s mysterious “‘first postulate’’ seems questionable; 
the presence of Example I (with its last sentence) im- 
mediately thereafter may lead to false notions. 

Chapter I, “Engineering dynamics and mechanical 
vibrations,” is the best. Chapter II, “. . . Tensor analysis 
of stationary networks and rotating electrical machinery,” 
(Kron’s theories) is the least convincing part. Chapter 
III, ‘‘Non-linearity in engineering,” considers: (1) methods 
of obtaining, or approximating, solutions to ordinary 
differential equations; (2) elliptic functions; (3) extensive 
bibliography. 

RICHARD S. BURINGTON 
Washington, D. C. 


The National Bureau of Standards 


The National Bureau of Standards, Washington, D. C., 
has announced two new mathematical tables which are 
now obtainable. They are the following: 
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Table of Fourier Coefficients (MT23). By ARNoLD N. 
LOWAN AND JACK LADERMAN. Reprinted from Journal of 
Mathematics and Physics, September, 1943, 11 pages. 
Obtainable from the National Bureau of Standards, 
Washington, D. C., for the sum of 25 cents. 

Coefficients for Numerical Differentiation with Central 
Differences (MT23). By HERBERT E. SALZER. Reprinted 
from Journal of Mathematics and Physics, September, 1943, 
21 pages. Obtainable from the National Bureau of Stand- 
ards, Washington, D. C., for the sum of 25 cents. 


Mathematics of Radio Communications 


By T. J. WANG. Pp. 371+ix, Figs. 199, 1422 cm. 
D. Van Nostrand Company, Inc., New York, 1943. 
Price $3.00. 


The vast expansion of the radio industry and its asso- 
ciated sister electronic industries, particularly in connection 
with the war, has stimulated a demand for a large number 
of technicians trained in the more practical aspects of 
laboratory work in electronics, but without the broad 
training of engineers. The principal obstacle in the way of 
advancement of these radio technicians is usually an 
inadequate training in mathematics occasioned by the lack 
of sufficient time to include the necessary formal courses 
in mathematics in their educational curricula. 

To remedy this situation, the book reviewed here 
sketches an outline for a course in the mathematics of 
most frequent use to workers in electronics laid out in a 
manner that may be given concurrently with classroom 
and laboratory work in electricity. The subject material 
has been well chosen and includes topics in arithmetic, 
algebra, analytic geometry, trigonometric functions, com- 
plex vectors as applied to circuit analysis, differential and 
integral calculus, and Fourier series. The examples and 
problems have been ingeneously chosen from many 
branches of electrical communication work, including such 
subjects as a.c. and d.c. circuit analysis, electron trajec- 
tories, vacuum-tube characteristics, analysis of oscilloscope 
figures, and the elements of modulation theory. This, 
together with a concise but well-written text, will no doubt 
contribute considerably toward making the book more 
interesting than the more usual mathematics texts to the 
group for whom it was written. 

Unfortunately, the lack of answers to any of the prob- 
lems will possibly inhibit the use of the book for self-study 
by a sizeable group of workers in the field, who no longer 
have the opportunity to attend classes, but who would 
have welcomed a volume of this type whose contents they 
might have mastered at home. 

LEsLiE L. FoLpy 
Columbia University, Division of War Research 


Basic Mathematics for War and Industry 


By Pau. H. Daus, JoHN M. GLEASON, AND WILLIAM 

M. WuyBurn. Pp. 277+ xi, Figs. 277, 14X22 cm. The 
Macmillan Company, New York, 1944. Price $2.00. 

On the home front, no less than on the fighting fronts, 

the present war has pointed out indisputably that the 
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basic mathematics ordinarily taught in high schools is not 
a dry academic subject limited in its usefulness to the 
torment of unfortunate students under the guise of 
“mental discipline,’ but that it is instead an active and 
useful tool for the solution of the most practical problems 
of modern warfare and industry. 

In view of the war-inspired need for rapid training, and 
review of large classes of our population in basic high 
school mathematics, the authors of the present book have 
tried to select the most useful portions of the usual syllabus 
in these subjects, and judiciously integrate them with 
their most important present military and industrial 
applications in a small volume which would prove equally 
valuable as a textbook in civilian and armed service 
training courses, as well as for home study. They have 
included in the volume material on arithmetic, algebra, 
plane and solid geometry, and plane and spherical trigo- 
nometry. There is a literal profusion of worked-out 
examples and problems (with answers supplied), the 
majority of which are drawn from machine shop work, 
air and sea navigation, physics, and similar war important 
applications. The entire emphasis is on mathematics as a 
tool, and the choice and organization of the material have 
been such as to expedite rapid familiarization with appli- 
cations, rather than appreciation of mathematics for its 
own sake. Thus, while the conventional form of presenta- 
tion of geometry as a series of propositions and construc- 
tions derived from axioms and postulates is retained, the 
logical aspects of the subject are completely submerged 
even to the extent of including considerably more axioms 
than are logically necessary. The text is lucidly and 
coherently written, and the typography is large and clear, 
making an attractive volume. 

This book is recommended to all who lacked the op- 
portunity to study high school mathematics or who were 
unfortunate enough to forget what they had once learned 
and who now find that such knowledge would be valuable 
to them in their present military or civilian occupation. 

Lesuiz L. Foipy 
Columbia University, Division of War Research 


Mr. Tompkins Explores the Atom 


By G. Gamow. Pp. 97+x, 17.5X22 cm. The Mac- 
millan Company, New York, 1944. Price $2.00. 
When discoveries have been made in physics by the 
application of new mathematical techniques to old prob- 
lems, there usually follow periods during which the full 
physical significance of the results is poorly apprehended 
by the majority of scientists. Gradually, while the new 
developments become more familiar by being discussed 
and applied in wider fields, their meaning grows more 
apparent, and simple verbal interpretations are found to 
supplement the original mathematical ones. For beginners, 
particularly, some mechanical model is usually the most 
acceptable, although once an initial period of professional 
apprenticeship has been served, it is the young men who 
go on to develop new ideas with confidence and enthusiasm. 
During this period of interpretation, it takes a man of 
unusual skill to create the necessary word pictures even 
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for a technical audience; and one of imagination, like 
Professor Gamow, to make a translation into the everyday 
language of the layman. In the latter case, an apparently 
fantastic story must inevitably follow. Perhaps this is why, 
in Mr. Tompkins Explores the Atom, Professor Gamow has 
proceeded with his hero’s education by giving him painless 
dreams which reveal the fundamental properties of the 
primary particles as though they were the natural in- 
gredients of a fairy tale. 

In three dreams, which occupy the first half of this little 
volume, Mr. Tompkins, or his wife, finds himself or herself 
on easy speaking terms with Maxwell’s demon, with the 
Electron tribe, and with the maker of Nuclei, from all of 
whom they receive intelligible explanations of the world 
of atoms. To drive his points home, the author adds in 
the second half of the book four lectures by ‘‘The Pro- 
fessor’ who illustrates by skillful analogy the behavior of 
atoms and their nuclei. A non-scientific reader will enjoy 
the fairy tales; a physicist, especially if he can recognize 
something of Rutherford’s features in the face of the 
woodcarver, will chuckle as he scans the other illustrations; 
and an advanced undergraduate student will find the 
dreams and the lectures useful as collateral reading during 
a first course in nuclear structure. Only a careless reader 
will suffer from philosophic insomnia by reason of the 
typographical error which introduces ‘‘fluctuations of 
destiny” on p. 10. 

If J. M. Barrie had been a scientist, M’Connachie, 
flying around lightheartedly on ore wing, dragging with 
him a reluctant author who would have preferred to be a 
family lawyer standing firm on his hearthrug among the 
harsh realities of office furniture—M’Connachie, instead 
of Mr. Tompkins, might have been the hero of a story 
such as this. 

Tuomas H. OsGoop 
Michigan State College 


Fundamentals of Optical Engineering 


By Donatp H. Jacoss. Pp. 487+x, Figs. 343, 
14214 cm. McGraw-Hill Book Company, Inc., New 
York, 1943. Price $5.00. 


This is a timely text suitable for a first course in the 
design of optical instruments. The wartime importance 
of this subject is brought home by the nature of the 
instruments discussed. This may cause much of the book 
to appear overspecialized in postwar times. However, the 
basic principles may be readily applied to broader fields 
by one who has mastered them. 

The book consists of four parts of which the first, 
consisting of 151 pages, presents an excellent summary of 
the principles of geometrical optics. The theory of stops 
is presented in a manner that makes this important subject 
easy to grasp by a beginner. This is followed by a discussion 
of representative instruments such as the various types of 
telescopes, gunsights, and rangefinders used in the military 
service. 

The third part, dealing with mechanical and electrical 
design, is a unique feature of this book. Many engineers 
differ in opinion as to whether optical design should precede 
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mechanical design or vice versa. The author observes that 
these are but two aspects of a single problem and should 
be worked out concurrently. He states that a primary 
purpose in writing the book is to show how mechanical 
and optical design are interrelated. Not only is this 
accomplished, but the student is given a good practical 
introduction into mechanical principles and machine-shop 
operations as they apply to optical instruments. It seems 
strange, however, that optical-shop procedures and ma- 
chines are not discussed while the mechanical ones are. 
Then again, optical tolerances are treated, but methods of 
inspection are not. No doubt the reason for these omissions 
is that their inclusion would make the book inconveniently 
large for a text on the fundamentals of design. 

The fourth and final section specifically deals with 
optical design. Ray-tracing methods and some short-cut 
formulae are presented. The theory is illustrated by 
application to the design of aplanatic objectives and 
eyepieces. The treatment is elementary but sound and can 
serve as an excellent introduction to advanced texts such 
as that by Conrady, whose notation has been adopted. 

The book is profusely illustrated by drawings and some 
halftones. There are numerous suggestions for further 
reading at the end of each chapter. There are enough 
problems to aid the student in testing his knowledge and 
in gaining mastery through application. This is an excellent 
textbook on a subject of great importance. 

JOSEPH VALASEK 
University of Minnesota 


Structure of Metals 


By CHARLES S. BARRETT. Pp. 567+xiii, Figs. 388, 
15X23 cm. McGraw-Hill Book Company, New York, 
1943. Price $6.00. 


Dr. Barrett has discussed in one book all the aspects of 
metals which have been studied by means of x-rays. These 
aspects include not only those covered in standard crystal- 
lographic treatises, but also those developed extensively 
during the past fifteen years. These latter include order- 
disorder phenomena, residual stresses, the changes oc- 
curring during cold working and the resulting structure 
and preferred orientations resulting therefrom, the pre- 
ferred orientations accompanying crystallization from the 
melt and accompanying recrystallization of deformed 
metals, the intermediate steps in age hardening and in 
phase transformation. The related subject of electron diffrac- 
tion is treated briefly. The discussion of this large range of 
topics is not confined exclusively to the findings of x-rays. 
Other pertinent experimental data are also introduced, as 
well as theoretical interpretations. 

It is believed this book is the first in English to describe 
clearly the method and use of the concept of reciprocal 
lattice. 

Dr. Barrett has had a long and varied experience in the 
application of x-rays to the solution of problems in the 
structure of metals. This rich background is reflected in 
the style with which his book is written. The general 
physics of each topic is discussed with a minimum of 
mathematics so that the book may be read fairly rapidly. 
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On the other hand, an exhaustive reference to the literature 
enables one to find readily more complete discussions of 
the details. 
C. ZENER 
Watertown Arsenal 


Magnetochemistry 


By Pierce W. SeELwoop. Pp. 287+ ix, Figs. 80, 
16X 23% cm. Interscience Publishers, Inc., New York, 
1943. Price $5.00. 


This book gives a useful survey of the field of magneto- 
chemistry, defined by the author as “the application of 
magnetic susceptibilities and of closely related quantities 
to the solution of chemical problems.” It covers portions 
of the field of magnetism different from the well-known 
books of Van Vleck and Stoner, and is a welcome supple- 
ment to them. The latter books are concerned primarily 
with the relations between the susceptibility and the 
structure of the atom or molecule, expressed in terms of 
spin and orbital quantum numbers and temperature. In 
Selwood’s book the author starts with these theoretical 
formulae and discusses the changes in structure that occur 
as one goes from one atom to the next in the periodic 
table or as one combines a given electropositive atom 
with a variety of electronegative ones, or as the groupings 
of organic radicals are changed. Magnetic properties are 
interpreted in terms of atomic structure as expressed by 
the quantum theory, and magnetic measurements are 
shown to provide a powerful tool for discriminating 
between various hypothetical structures. 

The field covered by this book is similar to that treated 
by Klemm in his Magnetochemie published in 1936. Since 
this date over 1000 papers on magnetochemistry have 
been published and the consideration given to them by 
Selwood makes his book a competent up-to-date review 
of the subject. 

The plan of the book may be stated briefly. A chapter 
at the beginning is concerned with methods of measure- 
ment and treatment of the data. A chapter at the end on 
“Applied magnetometric analysis’”’ mentions some appli- 
cations valuable chiefly for metallurgists. The other seven 
chapters forming the main body of the book are devoted 
to the diamagnetism and paramagnetism of atoms, mole- 
cules, and metals, and (briefly) to the ferromagnetism of 
metals and compounds. : 

An important part of the presentation of data is the 
free use of tables—showing, for example, the variation of 
susceptibility with position in the vertical columns of the 
periodic table—and of structural formulae for organic 
compounds. 

The book seems particularly adapted to the student 
who is planning to do research in magnetochemistry. It 
can be read with pleasure and profit by many others 
concerned with the interesting magnetic manifestations of 
atomic and mulecular structure, and the copious references 
to original sources enable the student to extend his reading 
with a minimum of effort. 

R. M. BozortH 
Bell Telephone Laboratories 
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Here and There 








The Establishment of a Division of High Polymer Physics 
in the American Physical Society 


Dr. Karl K. Darrow, secretary of the American Physical 
Society, recently announced that a Division of High 
Polymer Physics had been authorized by the council of 
the society, acting in response to a petition signed by 
thirty-one fellows and members. 

The field of this division is defined as the advancement 
and diffusion of knowledge of the physics of high polymeric 
materials such as rubbers, textiles, and plastics. Meetings 
of the division may be held (if the council so authorizes) 
in conjunction with or separately from meetings of the 
society. 

The council has appointed an “organizing committee” 
to take charge of the organization and the affairs of the 
division until by-laws are formulated and adopted by the 
council and a mechanism is set up for electing an executive 
committee. The organizing committee consists of F. G. 
Brickwedde, K. K. Darrow, G. B. Pegram, and A. E. 
Ruark, selected from the council; and R. B. Barnes, W. F. 
Busse, P. Debye, J. H. Dillon, W. J. Lyons, and L. A. 
Wood, selected from the members sponsoring the division. 
W. J. Lyons (Southern Regional Research Laboratory, 
2100 Robert E. Lee Boulevard, New Orleans 19, Louisiana) 
is serving as secretary of this committee. 

According to Article |X of the constitution of the society: 

1. The council may, upon petition by members of the 
society, form a division within the society charged with 
the advancement and diffusion of the knowledge of a 
specified subject or subjects in physics. 

2. Each division shall elect an executive committee, the 
chairman of which shall report its activities and needs to 
the council. 

3. Any division may be dissolved at the discretion of the 
council, 


National Electronics Conference 


Designed to fulfill the need for providing a common 
meeting place and discussion forum for technical personnel 
interested in applications in the broad field of electronics, 
the first annual National Electronics Conference will be 
held at the Medinah Club of Chicago, 505 North Michigan 
Avenue, Chicago 11, Illinois, on October 5-7, 1944. 
Emphasis will be placed on the engineering aspects of 
electronics in the conference which is described as a 
“national forum for electronic developments and their 
engineering applications.’ The conference is sponsored 
jointly by the Illinois Institute of Technology, North- 
western University, the Chicago Section of the Institute 
of Radio Engineers, and the Chicago Section of the 
American Institute of Electrical Engineers. Technical 
papers and discussion groups are planned to treat the 
communication, industrial, measurement, instrumentation, 
control, scientific, and medical applications of electronics. 

Those wishing to submit papers for presentation at the 
conference are invited to communicate with Professor A. 
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B. Bronwell, Chairman of the Program Committee, 
Technological Institute, Northwestern University, Evans- 
ton, Illinois, or Mr. Beverly Dudley, Secretary of the 
National Electronics Conference, 520 North Michigan 
Avenue, Chicago 11, Illinois. 


National Science Teachers Association 


Formation of the National Science Teachers Association, 
“to stimulate, improve, and coordinate science teaching 
at elementary, secondary, and collegiate levels of instruc- 
tion,” was announced by Dr. Philip G. Johnson, President 
pro tem of the new organization. It has been formed as 
the first step in a merger of two national science teachers’ 
organizations, the American Science Teachers Association 
and the American Council of “Science Teachers. It is 
expected to represent all teachers of science, and it will be 
affiliated with the American Association for the Advance- 
ment of Science and the National Education Association. 
Membership in the association will be open to all science 
teachers, and to anyone interested in science teaching. 
Provision is made for affiliation of other science teachers’ 
groups, and any such organization may affiliate. Head- 
quarters of the association, for the present, will be at 
Cornell University. 


The Charles L. Mayer Nature of Light Awards 


The National Science Fund announces two Charles L. 
Mayer awards for contributions submitted before January 
1, 1946 on the nature of light. A prize of $2000 will be 
awarded for an outstanding contribution to our basic 
understanding of the nature of light and other electro- 
magnetic phenomena which provides in terms intelligible 
to the community of scientists at large a unified under- 
standing of the two aspects of these phenomena which are 
at present jointly described by wave and by corpuscular 
theories. The second prize of $2000 will be awarded for an 
outstanding comprehensive contribution to a_ logical, 
consistent theory of the interaction of charged particles 
with an electromagnetic field including the interaction of 
particles moving with relatively high speeds. 

The first award is intended to encourage for the benefit 
of the non-specialist the interpretation of facts already 
known to the specialist. It is hoped that the second award 
will stimulate attack on one of the most fundamental 
unsolved problems in physics. In making these awards 
the National Science Fund will have the assistance of a 
special Advisory Committee consisting of Dr. E. U. 
Condon, Associate Director of the Research Laboratories, 
Westinghouse Electric & Manufacturing Company, East 
Pittsburgh, Pennsylvania; Dr. Karl K. Darrow, member 
of the technical staff of the Bell Telephone Laboratories, 
and Secretary of the American Physical Society; Dr. 
Robert A. Millikan, Chairman of the Executive Council, 
California Institute of Technology; and Dr. I. I. Rabi, 
Professor of Physics, Columbia University. 

Contributions may be submitted to the National Science 
Fund of the National Academy of Sciences, 2101 Constitu- 
tion Avenue, Washington 25, D. C. 
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News Notes 


Dr. John T. Tate, Dean of the College of Science, 
Literature, and the Arts of the University of Minnesota, 
who has been on leave since 1941, has asked to be relieved 
of administrative responsibilities. On completing his work 
with the National Defense Research Committee he will 
return to the university as Research Professor of Physics. 


Dr. Willis R. Whitney, honorary Vice President of the 
General Electric Company and first director of its research 
laboratory, has been elected an honorary member of the 
Electrochemical Society. The certificate of honorary 
membership will be conferred on him at the Milwaukee 
meeting. 


Hendrik Johannes Van Der Bijl, Research Physicist, 
Director of War Supply, Union of South Africa, Vice 
Chancellor, Pretoria University, has been elected a fellow 
of the Royal Society of London. 


The University of Wisconsin will confer the degree of 
doctor of science on Dr. Jesse T. Littleton, physicist, 
Assistant Director of Research and Development with 
the Corning Glass Works, Corning, New York. 


Earl Russell (Bertrand Russell), who has lived in the 
United States for some years, has been elected a fellow of 
Trinity College, Cambridge. He expects to return to 
England during the summer to resume his work in phi- 
losophy and mathematics. 


The Electrochemical Society 


The Eighty-Fifth Convention of the Electrochemical 
Society was held April 13-15 in Milwaukee. Dr. R. M. 
Burns, President of the Society, was in charge. The theme 
of the convention was “The electrochemist’s contribution 
to the war effort.’"’ The main technical session was on 
“powder metallurgy.” 

Honorary membership in the Society was bestowed on 
Paul J. Kruesi, President of Southern Ferro Alloys Com- 
pany, Chattanooga, Tennessee, and Willis R. Whitney of 
the General Electric Company, Schenectady, New York. 
The Young Author’s Cash Prize and Book Prize go to 
Walter G. Berl of Carnegie Institute of Technology. 


The Physics Club of Philadelphia 


At a recent meeting of the Physics Club of Philadelphia 
held at the University of Pennsylvania, Dr. Martin A. 
Pomerantz of the Bartol Research Foundation gave an 
address on ‘‘The nature of cosmic-ray particles.” 


Necrology 


Dr. Margaret E. Maltby, Associate Professor of Physics, 
retired, of Barnard College, Columbia University, died on 
May 3 at the age of eighty-three years. 


Dr. Lionel Robert Wilberforce, Professor of Physics at 
the University of Liverpool from 1900 until his retirement 
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in 1935 with the title emeritus, died on April 1 at the age 
of eighty-two years. 


Erratum 


Through a regrettable accident, part of the heading in 
the advertisement of Ballantine Laboratories on page ii of 
the May issue was dropped out in the printing. It should 
have read 0.00002 to 10,000 volts. 





New Booklets 








The Harshaw Chemical Company has released a book- 
let entitled Synthetic Optical Crystals which describes 
the development of the large scale growing process and 
the production of synthetic crystals. A limited number of 
booklets are offered free upon-request by The Harshaw 
Chemical Company, 1945 East 97 Street, Cleveland, Ohio. 


Nickel Steel Topics, Vol. 13, No. 1, March, 1944, is an 
8-page booklet which is published in the interest of pro- 
ducers and users of nickel alloy steel alloys. Requests for 
copies should be addressed to The International Nickel 
Company, Inc., 67 Wall Street, New York 5, New York, 
attention of H. S. Lewis, Editor. 


Chemical Digest, Vol. 10, No. 1, Spring, 1944, contains an 
article entitled ‘‘Chemical products for postwar needs” as 
well as other items concerning chemical news. Copies of 
this 4-page publication are available upon request from 
Foster D. Snell, Inc., 305 Washington Street, Brooklyn 1, 
New York, attention F. E. Parry. 


Wheelco Comments, Vol. 3, No. 3, January-February, 
1944, is a 4-page booklet which presents articles.on the 
problems of industrial measurement and control. It is 
published by Wheelco Instruments Company, Wheelco 
Building, Harrison and Peoria Streets, Chicago 7, Illinois, 
and available on request 


A new folder which describes the new model Spencer 
refractometer is available and will be supplied by the 
Spencer Lens Company, 19 Doat Street, Buffalo 11, 
New York. 


The new Cannon Electric Laboratory and Switchboard 
Connector Bulletin contains information, photographs, 
drawings, and data on this line of specialty fittings. 
Catalogued are surface and sub-mounting plugs and 
receptacles, straight cord plugs and receptacles, switching 
plugs, and experimental switchboards. This 12-page 
bulletin will be sent upon request from Cannon Electric 
Development Company, 3209 Humboldt Street, Los 
Angeles 31, California. 


The Ohmite News, April, 1944, contained a résumé of the 
career of Heinrich Friedrich Emil Lenz originator of Lenz’s 
law. Described also is the Ohmite rheostat and its appli- 
cations. Ohmite Manufacturing Company, 4835-41 Flour- 
noy Street, Chicago, Illinois is the publisher (2 pages). 
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Transmission Line Theory Applied to Wave Guides and Cavity Resonators. I* 


Davip MIDDLETON AND RONOLD KING 


Cruft Laboratory and the Research Laboratory of Physics, Harvard University, Cambridge, Massachusetts 


In Part I circuits consisting of hollow metal wave guides 
of arbitrary cross séction are analyzed in terms of scalar 
and vector potential and stream functions in a form 
representing a generalization of the treatment of the 
conventional lines given by R. King. Completely hyper- 
bolic solutions analogous to those given in the usual 
theory are obtained for the E and H modes; these yield 
general relations for the stream and potential distributions. 
Terminal, input, and characteristic impedances are defined 
for the wave guide or cavity resonator. The propagation 
constant is discussed, and it is demonstrated that the 
attenuation (constant) due to the imperfect metal walls 
of the pipe and the dielectrics in the guide are additive if 
the materials have low loss. The characteristic impedances 


(Received February 1, 1944) 








for the two types of mode are considered, and general 
formulae for the practical conditions of operation are 
included. The analogy between, this theory and that of 
the conventional line is discussed, and it is pointed out 
that the physical interpretations of results formally 
identical in the two cases are equivalent for a given mode 
and a particular cross section of pipe. The analogy breaks 
down, however, when one attempts to treat different 
modes and guides of different transverse dimensions at 
the same time, for example, in the problem of impedance 
matching. One advantage of the theory is that experi- 
mentally meaningful results may be obtained for general 
terminations which ordinarily cannot be handled ana- 
lytically. 





INTRODUCTION 


HE general analysis of the conventional 
open-wire and coaxial transmission lines, 
closely spaced compared with a wave-length, is 
well known, and has been expressed in many 
different forms. It has also been shown that the 
fundamental differential equations relating the 
total current J and the potential difference V for 
the conventional line are formally analogous to 
the relation governing the behavior of metal 
pipes and coaxial tubes whose cross sections are 
sufficiently large to allow higher modes in the 
distribution of charge and current. This formal 
analogy is obtained by replacing the current J 
and the voltage difference V of the conventional 
line with a scalar (or vector) stream and po- 
tential function, respectively, from which it is 
possible to determine the fields and the charge 
and current densities. 
_ Perhaps the most convenient general solution 
of the transmission line equations is the com- 


*A considerable part of Part I is a modification of 
course lectures given by the second-named writer at Har- 
vard University during the last three years. 
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pletely hyperbolic one, given by R. King,’ in 
which the various terminal impedances are 
expressed in terms of hyperbolic functions. The 
main advantage of this form is that the over-all 
effects of the attenuation and phase shift of the 
terminations may be added, respectively, and 
directly, to the over-all attenuation as and phase 
shift Bs of the line itself, and that the solutions 
are quite general, and in the practical cases yield 
simple results for the determination of such 
quantities as input impedance, standing wave 
ratios, transfer of power, etc. 

It is our purpose in this paper, which at least 
in Part I is to be regarded as a companion piece 
to reference 1, to show in detail that the analysis 
of transmission circuits constructed of hollow 
and coaxial conductors may be given in the 
general hyperbolic form, mentioned above, of 
the ordinary theory, and to apply that theory 
in a number of practical cases. Specifically we 
will define and discuss terminal and input 
impedances of a wave guide and cavity resonator 





'R. King, J. App. Phys. 14, 577 (1943). 
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of arbitrary cross section, unrestricted termina- 
tions, and driven at some arbitrary point along 
the tube. A brief discussion of the propagation 
constant and characteristic impedances will be 
included, and finally in Part II we will consider 
a number of experimental examples, and a 
theoretical treatment of a closed, cylindrical pipe 
containing several imperfect dielectric media at 
one end. The great advantage of the analysis 
lies in the fact that definite and meaningful 
experimental results may be obtained, for any 
termination, whereas a mathematical treatment 
of the problem would, in most cases, be quite 
impractical if not impossible. 

Since, as will be shown, there exists a formal 
analogy between our results and those of the 
conventional theory, the formulae obtained in 
one case apply equally well, formally, in the 
other, so that it is a comparatively simple matter 
to obtain results in either system, once the 
necessary electrical and geometrical constants 
of the particular transmission circuit under 
examination have been measured. For this reason 
we will attempt to follow the notation of a 
recent article on transmission lines,! and will on 
occasion quote formulae rather freely from it, to 
speed the use of our derivations in practical 
problems. It must be emphasized, however, that 
the analogy between the two cases is not com- 
plete, and hence care must be taken in the 
interpretation of the results obtained in one case 
when applied in the other. We will consider the 
nature of the analogy in a more detailed manner 
at the end of Part I. 


PART I. DEVELOPMENT OF THE THEORY 


In deriving expressions for the stream and 
potential distributions, from which we will make 
our definitions of impedance, we proceed first 
from Maxwell’s equations for media free of 
sources and sinks: 

Vv XE= —iwuH = —ZnH ’ 
Zm=twu, (ohms/m), (1) 


Y,=iwe*, (ohm—m)~!. 
VXH=iwe*E= VE; (2) 


e* =€—10/w. 
Vv -H=0, 
v-E=0. 
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These relations are the fundamental field 
equations for the steady state, or forced oscilla- 
tion phenomena such as we are dealing with 
here, for which the time variation of the fields is 
of the form e“‘, where w is the angular frequency 
of the source. No restriction is placed on the 
various media other than that they be simple, 
homogeneous, and isotropic, the usual conditions 
in practice. Then the (absolute) permittivity e, 
permeability y, conductivity o, measured in 
farads/m, henrys/m, and mhos/m, respectively, 
are constants independent of the applied and 
induced fields. Because of the solenoidality of 
the E and H fields, vide Eqs. (3), (4), vector 
stream and potential functions may be defined, 
along with a pair of scalar functions which are 
chosen in such a way that the vector quantities 
are uniquely specified. That this may be done 
follows at once from the definition of a uniquely 
determined vector,? which requires that its curl 
and divergence be known, and that it vanish at 
infinity as 1/r?. As we shall see presently, the 
pair of functions derived from V -E=0 gives rise 
to transverse electric (TE),* or H-type modes, 
and those from V-H=0 yield the transverse 
magnetic (7.M)* or E-type modes, where the E 
and H types are distinguished here from one 
another by the superscripts g= 1, 2, respectively. 

In order to have available results to which we 
shall wish to refer later let us proceed in a 
somewhat more detailed fashion than we would 
ordinarily. Considering the HT modes first, we 
define a vector potential A® in volts from (4) 
as follows: 


V XA? = —E), (5) 


since the operator ¥V-V xX always yields zero. 
From (2) this gives 


VX(H®+ YA) =0, (6) 


from which a scalar stream function ¥® (am- 
peres) may also be defined, so that 


vy =—(H+ VA), (7) 


whereupon we obtain 
H® =— Y,A®—yv), (8) 


2M. Mason and W. Weaver, The Electromagnetic Field 
(University of Chicago Press, 1929), theorem IV, page 357. 
3S. A. Schelkunoff, Electromagnetic Waves (D. Van 
Nostrand Company, Inc., New York, 1943), pp. 153, 154. 
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Now 

VXVXA®=-VXE® =Z,H® 
=—Z,(VYA®2+yV¥). (9) 

In order to be able to separate the variables and 


to determine A® uniquely, we relate ¥® and 
A® by the following equation 


Vv -A® = —Z,,v), (10) 


From (9) and (10) and the vector identity 
V-VA=VV-A—YVxXVY XA we obtain finally the 
familiar vector and scalar wave equations for 
the transverse electric modes: 


V-VA%+hA™ =0, (11) 
VO+kRYO=0, g=1, 2, (12) 

where 
k?= —Z,, Y.=w*pe*. (13) 


In a similar way, because of the vanishing of the 
divergence of H, we may define for the E modes 
a vector stream function A® (amperes) and a 
scalar potential ¥® (volts). The uniqueness 
condition (10) now becomes 
VAM=—-—YV,V", (14) 
and Eqs. (11) and (12) are seen to hold equally 
well in this case. The fields are given by 
E® = —Z,AY-— VV, 
HM=y7xXxA". 


(15) 
(16) 


The vectors AM, g=1, 2, may in general have 
three components ; however, this is a considerable 
and unnecessary complication, as the vector 
operation ¥V-¥V does not admit the solution of 
Eq. (11) by the usual method of separation of 
variables for non-Cartesian coordinate systems. 
This difficulty may be avoided if we specify that 
A“ has only one component. No generality is lost ; 
the fields of a particular mode are still completely 
determined. The most general expression of the 
fields requires a superposition of both the E (¢=1) 
and H (q=2) types of oscillation.‘ This require- 
ment is true, regardless of the number of compo- 
nents of A™. For axially symmetric systems, 
such as the ones considered in this paper, the 
natural choice of the components of A® is 
obviously that in the z direction, i.e., 


A‘ =i;A,, g=1, 2, 


*E. Ledinneg, Ann. d. Physik 41, 557 (1942). 


(17) 
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where i; is a unit vector in the 2 direction. With 
this in mind, and separating the variables, we 
may write (11) and (12) as 
2 
Vu Ag’? +(y@)2A,@ =0, 


07A,“@ 
——(K,)?A,@ =(), 


(18) 





g=1,2, (19) 


2 


Oz 


and similar relations for V‘®. Here v., is the 
Laplacian in the curvelinear coordinates u, v of 
the tube cross section. The quantity y™ is the 
“shape factor,”’ and is (in general) complex if 
the walls of the wave guide are not perfectly 
conducting; otherwise 7 is real. From the 
separation of variables we obtain the following 
relation between K,@ and y™: 


(Ky)?= (7)? k= [ay +18, F, 


and as we shall see shortly A,“ may be identified 
as the propagation constant of the pipe. If we 
use (10), (14), and (17), and similar relations 
for the 7M modes we find at once the familiar 


(20) 


equations 
0A,‘ . 
a —=—( yn, (21) 
dz Zam 
anc awa y, 
ae cxwopase/(), en 


where the upper quantity in braces applies for 
g=1 and the lower for g=2. These equations 
are seen to be formally identical with those of 
the conventional line,® provided we define char- 
acteristic impedances for the E and IH modes by 


2.2K, /¥Y., (23) 
Lem=Zm K,®, (24) 
and provided we identify the potential and 
stream functions of a particular mode with the 


corresponding. potential difference V and the 
total current J of the usual theory. (See Fig. 1.) 











YY” 


ee < e > = 
» @ ) @) 
AO | A Wt 4 


2 4 s-z 





Z=s 


Fic. 1. Circuit diagram of a wave guide and 
cavity resonator. 


5 See reference 3, Chapter VI, Section 11. 
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SOLUTIONS OF THE AXIAL EQUATIONS AND 
DEFINITION OF THE TERMINAL 
AND INPUT IMPEDANCES 


II 


Typical solutions of the axial equations (19) 
are* 


cosh whee, 
sinh K ©)z 


Ay (u, veto! -|cu| 


{sinh Kz 
{@ 


+D | [rec sdeis (25) 


lcosh Kz 
where the upper and lower lines apply when 
qg=1, 2, respectively. The functions T“ (u, v) are 
solutions of the cross-sectional equations (18) 
and in the next section will be considered briefly 
for two special cases. Relations similar to (25) 
may be obtained for ¥™ at once from (21). The 
constants C; and D,™ will now be determined 
from the boundary conditions at z=0,s. Con- 
sidering the JZ modes first we find that the 
boundary conditions on (19) are, at z=0, 


A, (0) = Ao® (volts) (26) 


for which A»® is not specifically known in 
general, and is a function of the source, detector, 
and any other apparatus present in the pipe. 
lf A. (volts) represents the potential distribu- 
tion of the source for a given mode we have 


A, =A,)(u,v) =B,T (u,v), (27) 
and then, from Fig. 1, we see that 
Ay SA, —V9Z_, (volts), z=0. (28) 


This defines the terminal impedance Z,)®, which, 
like Ap and W,®, depends on the source, its 
coupling to the guide, etc. From (10), or what 
is the same thing in view of (17), from (21), we 
obtain for (28) 

Zo OA 1°?) 


A (2) = 4 (2) +4— eunnncinienia 29) 
, | a Oz z=0 


where we specify in general that 
Ao) =ByT®)(u, v) ; (30) 


here By® is a constant depending on the ampli- 
I 


* We have omitted the subscripts g on A, ¥, K, etc., 
whenever there is no ambiguity in so doing. 
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tude of the source, and Eq. (30) implies that our 
termination at z=0 is sufficiently distant, a 
wave-length or so in practice, from the exciting 
probe (antenna) so that the transverse distribu- 
tion of the vector and stream functions is that 
of the natural mode excited, i.e., is the solution 
of (18). We will not attempt to consider the 
near-zone development of the fields, which is at 
best an extremely laborious problem, for the 
far-zone theory treated here is in most cases of 
much greater practical importance. 

For the arbitrary termination at z=s we define 
the impedance by 


Z,;° =A,°)(s) /W,(s) 
Ai! 
= —ZmnAx(8) [| (ohms). (31) 
2 lz=s 





By satisfying the appropriate conditions for 
E® and H® at the boundary z=s, which in turn 
requires that we can determine the electro- 
magnetic field distributions for the termination, 
we obtain an explicit analytical expression for 
Z.™. Usually this is not possible; but it will be 
demonstrated in Part II for the special case of a 
resonator or guide terminated in several trans- 
verse sections of dielectric material. 

In an analogous fashion we find for the E 
modes that 


WoO = VY —Z,Ao™ (volts), (g=0). (32) 
This defines the terminal impedance Z), similar 
to (28). From (14) it follows immediately that 

dA, 
#0 — ZAM = — Ze - £2,21/Y.. (33) 


2 |e=0 
Also we have 


Ay (0) =Ap =BoT(u, v), 
Vv. (u, v) = BOT (u, v). (34) 


and 


At z=s the terminal impedance for the 7M 
modes is defined as 


Zi) =W,0(s)/A1M(s) 
0A,°) 


-s 


/ Av) (ohms). (35) 
With these definitions, (28), (31), (32), (35), as 
one would expect, our problem becomes formally 
analogous to that of the conventional line, 
schematically represented in Fig. 1, with the A’s 
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and ¥’s, depending on the mode excited, replaced 
by the appropriate I’s and V’s of the usual 
theory. The final results are then expressible in 
the completely hyperbolic form,' and we have 
for the H modes 


W,® (amp.) 
A, (volts) 


A, (u, v) sinh 69° 
~ sinh (Ks+6)2)+0,) 
| 1/Z-m sinh (K®w+é,) 
| cosh (K@w+0,)) 











; (36) 


where the upper line in the braces applies for 
v,, and the lower for A,®. For the E modes 
we obtain similarly 


A,“ (amp.) ¥.(u, v) sinh A 
sa (volts) | ” sinh (KYs+69 +0,) 
1/Z¢ cosh (K“w+6,“) | ; 

| sinh (K®w+0,() I" (37) 








Here w is the distance from the load or termina- 
tion Z,™ to the point z where A, and ¥, 
are computed. Thus 


w=S—Z. 


(38) 


The time variation may be simply introduced by 
multiplying (36) and (37) through by e“'. The 
complex quantities 0), 0, are defined in terms 
of the p and @ functions! by 


Oo = po? +14, 6,°% = p,(*%) +716, “, 


and the terminal impedances are related to p 
and @ in the following manner: 


(39) 


ce 


Zy\0 =| coth @)‘@ ; 


rr (40) 


Z,0 =} coth 6,°®. 





cm 


If p, and #, are known, for example, from 
experimental data, we may obtain the terminal 
resistance and reactance at z=s from 


sinh 2p,” —@.{%&} sin 26,‘ 


cosh 2p,‘ —cos 2®,‘® 





Ri =Rel] | (41) 


and 
X,Y = —R.i% 
>" 24,‘ +¢.{%,} sinh ee) 


cosh 2p,‘% —cos 2,‘ 





(42) 
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Here the upper line in the subscripts {4,} is to be 
taken for the E modes and the lower line for the 
H. Reciprocal relations for p,“ and ®, in terms 
of R,® and X, may be found at once from (41) 
and (42). We have then 


ps‘? = tanh"! 
2Re\ R,“? -— sits X,&@ 
tim} ( bet in} ) | (43) 
RO?FXOFRE | (1t Gel} 9) 





and * 
&,(9 =} tan— 
—2Ref in} (Xe +Ge a} Ro 
: } ( Gen} ) | (44) 
R,O24X,92?-—R2f 5} (14+ef%} 2) 


&,=2/2 for a metal piston. 





It is often convenient to define the dimensionless 
ratios 7; and x; by writing 








ZO /Zek om} =r +1x1., (45) 
where, following King,* we have 
R.© — del in} Xe 
71.6? =——_ _— = 
Rel im} (1+ Gc! in} *) 
sinh 2,‘” 
=- —_———_———., (46) 
cosh 2p,‘?) —cos 2,” 
and 
(q) e (@) 
X,(ou-* +¢.{%}R, 
Re} in} (1 +c in} ?) 
—sin 26,‘® 
= (47) 





cosh 2p,‘ —cos 26, (a) 


In all but a few critical cases the terms involving 
¢-{%} are negligible for good guides, and so we 


have 
Pr is6? = R,©/R-. ; ‘) ? 


(48) 
° X15? = X,, Rett} . 


The definitions (40) may be extended immedi- 
ately to include the input impedances in the 
following way: 


ce 


Z:,°%= coth (A,“+72F,), (49) 





cm 





* Here and in reference 1 the following convention is used: 
, isin quandrant I if the signs of the argument of the arc- 
tangent in (44) are +/+; in quadrant II if the signs are 
+/—; in quadrant III if the signs are —/+; in quadrant 
IV if the signs are +/—. 

6 Reference 1, formulas (19a) and (19b). 
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where 
AM =a%stp.; FO =B%s+,, (50) 


Then we may obtain the input resistance and 
reactance for the E and H modes from (41) and 
(42) simply by replacing the subscript s by in 
and p.“”, ©,” by A, and F,™, respectively, 
and similarly for the reciprocal relations (43) 
and (44). 

The expressions developed by King,! formulas 
(21)—(71), including Table I and Figs. 6, 7, 12 of 
that paper, may be applied equally well in the 
corresponding analysis of wave guides given 
here, so we will not attempt to carry the details 
of the discussion of impedances and admittances 
further. 


POTENTIAL AND STREAM DISTRIBUTIONS 
III 


Before examining the axial behavior of A,“ (z) 
and W,(z), let us consider the transverse 
equations (18). Particular solutions may be 
obtained at once for the coordinate systems for 
which the variables (u,v) are separable.’ For 
the two particular cases of a rectangular and 
circular cross section, respectively, we have the 
following possible solutions, from which we may 
construct the proper one to fit the transverse 
boundary conditions below : 

Sin d(x} (cos b,°%y 
To(u, 0) =| | ; ¥ 
COS dm'?x} (sin b,%y 
_m=0, 1, 2, 3,--- 
n=Q, 1, 2, 3, --- 
(@) 


(dm)? + (bn°)? = (ym.n)?; OSxSa; OSy=d; 


(S1) 


and for the circular pipe of radius a 





= (q) . | COS n@é 
1 (u,v) =Jn(Ynpr) ? ’ 
sin 20 
n=(Q, 1,2, 3,--: 
(52) 
p=1, 2, 3, --- 


OSrse; 050222. 


Here J, is a Bessel function of the first kind and 
order n. The functions listed in (51), (52) are 
sufficient for our discussion and for the analysis 


7J. A. Stratton, Electromagnetic Theory (McGraw-Hill 
Book Company, Inc., New York, 1941), p. 349. 
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of the problem treated in Part II. [For a more 
complete list of solutions to (18) in the circular 
case, the reader is referred to any standard work 
on differential equations. ] 

The boundary conditions to be satisfied at the 
inner walls of the hollow pipe are the familiar 
ones® requiring, for homogeneous and isotropic 
media, that the tangential component of the E 
fields and normal component of the H fields 
vanish, i.e., 


nXE,=0; n-H,=0; D, =«Ei; B, =H. (53) 


If » is the surface charge density (coulombs/m) 
and I the surface current density (amperes/m), 
then for perfectly conducting walls 7 and I are 
determined from 


n-Di=n; nXHi=I. (54) 


The volume density of current is J=cE (amp./ 
m*), and exists only in the dielectric, provided, of 
course, that Gaiciectrie#O. The subscripts (1) on 
(53) and (54) denote the interior of the hollow 
pipe. If the walls are not perfectly conducting, 
however, the fields do not vanish at the inner 
surfaces of the tube or within region (2), and 
(53) and (54) become 


n X (E,—E,) =0; n-(B:—B.)=0, (55) 
n:(D,—D:)=n; nX(H:i—H2)=0. (56) 
Using the asymmetrical form of Green’s theorem 
applied to a bounded surface,®!° and with the 


aid of (18), we obtain the following relation for 
the ‘‘shape factor” y of the hollow pipe: 


f ftv. yl (o* ‘Tin. T |dS 
; 8 


- f roa -VT)ds 


(y0)?=— , (87) 


[ [reresas 


8 





where S is the cross section of the wave guide 
and s its periphery. For perfectly conducting 
walls n- V7 or T™* vanish, as a consequence 
of (53), and (7)? is seen to be real. However, 
8 Reference 7, pp. 34-38. 
® Reference 3, Chapter 10, Section 4. 
10 R. I. Sarbacher and W. A. Edson, Hyper and Ultrahigh 


Frequency Engineering (John Wiley & Sons, Inc., New 
York, 1943), Chapter VI, Sec. 19, Chapter VII, Sec. 14. 
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if the walls are not perfectly conducting, the line 
integral in (57) does not become zero, because 
of (55), and (y“)? is then complex. 

The general problem of hollow pipe trans- 
mission may be considered in two parts, one, in 
which the source is closely coupled to the wave 
guide, or (in some cases) to a cavity resonator, 
and the other, in which the source is loosely 
coupled, approaching the free oscillation condi- 
tion. For the transmission of power we are 
primarily interested in close coupling and in 
having the generator at one end of the system 
over which energy is being propagated. The 
solution of this problem has already been ob- 
tained in the previous section, and the final 
equations of the potential and stream distribu- 
tions are given by (36) and (37). Since we have 
strong coupling we cannot speak of “point” 
sources, in the conventional transmission line 
case, or correspondingly, of ‘plane’ sources for 
the hollow pipe, but rather we must specify a 
volume source distribution, whose exact nature 
it is rarely, if ever, within our power to describe 
accurately. For this reason we have not at- 
tempted to analyze the fields, charge, current, 
and potential distributions in the immediate 
neighborhood of the source, but have assumed 
that at a certain distance from it, and any other 
apparatus present, detectors, etc., of no more 
than about a wave-length in the closely coupled 
cases, and somewhat less for loose coupling, the 
transverse field pattern of the particular mode 
excited is known, and corresponds to that of a 
natural mode of the tube. This procedure 
eliminates the intricate and approximate analysis 
which attempts to deal with source geometries, 
and is reasonably accurate even if the generator 
does not have the same field pattern in its near 
zone as it does at a distance of about a wave- 
length. It is important to note (and it will be 
considered more fully at the end of Part 1) that 
for the wave guide or resonator theory the 
potential and stream distributions are three- 
dimensional, comprising a cross-sectional as well 
as longitudinal variation, whereas the current 
and _potential difference are one-dimensional in 
the conventional transmission line analysis. 

It is possible to generalize the solution obtained 
in the previous section to account for a volume 
source distribution at any region along the tube. 
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This is not the problem of greatest practical 
interest, however, since generally for power 
transfer it is more usual to have the source at 
one end, feeding a load at the other, a problem 
already successfully treated. 

For the case of loose coupling we can consider 
the’ generator to be a plane source, whose field 
distribution is known over a cross section of 
infinitesimal width. Such a distribution may be 
physically realized by a series of antennae 
suitably arranged. A single antenna will also 
serve, but we cannot assume that the relations 
for the fields, charges, etc. will hold to within 
the same distance of the source as for the more 
elaborate distribution. 

The loosely coupled case may then be easily 
generalized for a source at any point along the 
tube." If the guide (or resonator) is driven by a 
loosely coupled oscillator, which sets up a cosine- 
symmetrical field with respect to the planes 
z=x,"'8 along the pipe, then it is sufficient to 
replace 6) in the numerator of (36) and (37) by 
6° +Kx, and to write s—x, for s, with s the 
new over-all length. 

We then have 





sinh (K“s+6)2+0,) 


Vv, Ac.2,(t, v) sinh (K®x,+ 69) 
\ ~ 


1/Z.m sinh (K®w+é6,) 
, (58) 
cosh (Kw-+ 6, ) 
for the H modes, and for the E modes 
[Ax | We, 20(t, v) sinh (Kx,+69) 
Ly, | ys sinh (Ks+6) +4,) 
1/Z.. cosh (Kw+,) 
(59) 


sinh (KMw+@,“) 


If the loosely coupled oscillator sets up a sine- 
symmetrical field, so that it is equivalent to a 
“double layer” plane distribution of sources, 
symmetrically placed with respect to x,, we have 
results similar to (58) and (59), except that 


sinh (K‘x,+ 69‘) 


‘For the conventional line see R. King, Proc. I. R. E. 
27, 715 (1939). 

2 See also R. King, Proc. I. R. E. 28, 84-87 (1940). 

1 R. King, Rev. Sci. Inst. 10, 325 (1939). 
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CS 





is replaced by 
cosh ere 
(1) 


Pe me and ve A by Fe . and We, Zs. 
For an asymmetrical oscillator a sum of dis- 


; . (2) — (2) 
tributions using Ae,z, and Ae,z, for the H modes 


(1) — (1) 
and WVe,r; and We,z, for the E modes may be em- 
ployed. Cosine-symmetrical distributions may be 
obtained by placing the antenna, in the form of 





a circular loop at a distance of n half-wave- 
lengths in the guide from one end of the tube. 
When a straight wire antenna is located an 
odd number of quarter-wave-lengths from one 
end, sine-symmetrical fields will result.14-!7 

In practice we are usually concerned with the 
magnitude of A; and ¥,, which become for 
the case of a source at one end, from (36) 


and (37): 





ey" |=14 eo 
[4 


|| Ay 


sinh? po’? +sin? Bo‘ j 
sinh? (a!s-+ po® +p.) +sin® oper 
1/|Zem| [sinh? (aw p,) +sin? (B@w+,) ]! 
| | [sinh? (aw p,“) +cos? (8Bw+,) J} 


. 


sinh? po’? +sin? By? 





The condition for resonance (maximum values 
of the amplitudes) is in all cases 


BOs, +-6)°9 +6, =n, (62) 
provided 


[(a /B™) sinh (a s+ po + p,“?) 21. (63) 


The width As, of a resonance curve between 
the lengths for which the amplitudes |¥|, |A|, 
etc., have 2~'/? of the maximum value may be 
obtained from (60) and (61). It is 


As, =2(a's, + po +p.) /B™, (64) 


provided 
(As,(® /2)?<(s,°)?. (65) 


THE PROPAGATION CONSTANT 
IV 


The propagation constant K, has been given 
by Eq. (20), where k?=wye* [vide (2) ]. Let us 
examine the shape-factor (y), which by (57) 
is in general complex. As we have seen from 
(55) and (57) this complexity is the result of 
the finite conductivity of the metal walls of the 
pipe. Omitting the superscript (q) for brevity, 
since these results hold for either type of oscil- 
lation, we may write 


y= (yret+1yi)?=77(1+1h.), (66) 
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sti Pe (a s+ po +p.) +sin? Gee 
| 1/|Zem | (sinh? (a®w+ p.“) +cos? (BOw+e) ]Aj 
[sinh? (a? w+ p.) +sin? (BPw+eh, }} 


' (61) 





where : ‘ 2 
=y2—y72; he=2yryi/F?. (67) 


Then we have 
K,?=(72—Ba®) +i(hey?+haBa®), (68) 


Ba? =wpaea; ha=ca/wea, (69) 


and 


where the subscript d denotes the dielectric and 
the subscript c refers to the conductor. If we 
make the further definition 


vr =7*/Ba° (70) 


then (68) becomes 


hath.v* 
K,?= —8a(1—v? »(1-|~—=}). (71) 
1—p? 


We have also the following relations between 
the phase and group velocities v, and vg, and 


“LL. J. Chu and W. L. Barrow, Proc. I. R. E. 26 (De- 
cember, 1938). 

15 G. C. Southworth, Proc. I. R. E. 25, 807 (1937). 

16 Reference 10, Chapter VIII. 

17 J. G. Brainerd, G. Koehler, H. J. Reich, L. F. Wood- 
ruff, Ultra-High Frequency Techniques (D. Van Nostrand 
Company, Inc., New York, 1942), seventh printing. This 
book contains an extensive bibliography in Chapter 16. 
See also a Reference Guide to Ultra High Frequencies, ¥. 
Kelsey (Zenith Radio Corporation, Chicago, 1942); E. l 
Condon, Rev. Mod. Phys. 14, 341 (1942). 
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the wave-length in the guide \,: 
Vg=w/B,; Va =dw, OB,; 


catia (72) 
g— 47/Aqg; 


UVg0G =Ca’. 

Here cy (=paea)~') is the velocity of propagation 
of plane waves in medium d and 8, is the phase 
shift, obtained from (20) and (71) and given 
explicitly below. We have three cases to distin- 
yuish. 


Case I. Transmission, v?<1: From (71) we 
have 
K,=a,+i8, =Ba(1 —v?)*Lg(h) +if(h) ], 


(73) 
h=(hath.v?)/(1—v?), 


where f(h) and g(h) are functions tabulated 
originally by G. W. Pierce,'* and are defined as 


f(h) ={C(1+h?)!4+17/2}! 


=1+h?/8—5h*/128+---, (74) 
g(h) ={((1+A?)!—1]/2}} 
= (h/2)(1—h?/8+9h*/128—---). (75) 


From Eq. (73) we obtain the general expression 
for the attenuation and phase constants 


8B, =Ba(1—v?)' f(h). 


For perfect conductors and dielectrics h.=ha 
=0,. and consequently g(h)=0; f(k)=1. With 
these values it is a simple matter to determine 
a,, By, v, Us, Ue from (76), (72), (70), (67). By 
far the most important case is that for low 
dielectric losses and good conductors. This re- 
quires that 


at, = By(1 — v?)!g(h); (76) 


h?/8&<«1, (77) 

as can be seen at once from (74). Then we have 
Bav*h, Baha 

Wy = Age F Od =; Torre (78) 


2(1—v%)#  2(1—»2)!” 


and it follows immediately that (a,/8,)’<1, 
which is the usual condition for good transmis- 
sion circuits. It is quite valid to consider the 
attenuation as a linear combination (78) of the 
separate attenuations due to the dielectric and 
tube walls, provided the condition (77) is ful- 
filled. This linearity is commonly assumed in 


8G. W. Pierce, Proc. Am. Acad. Arts and Sci. 57, 
175-191 (April, 1922). 
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practice but does not seem to have been clearly 
demonstrated previously. For poor dielectrics 
and metal walls higher powers of h, cf. (74), are no 
longer negligible, and the linear relation (78) is 
not then correct. In this case we must use (74) 
and (75), with the aid of Pierce’s tables.'’ 

For good dielectrics and conductors we have 


QAygd = (Ma €a)*Oa 2(1 —vp*)}; 
By =Ba(1—v?)}; 


Vs =Ca(1—v?)-3; 


p (79) 
v=%¥r"/ Ba"; 


Vg =Ca(1—v?)!. 


(80) 


Case II. Cut-Off, v=1: In this instance (68) 
reduces, with the aid of (70), ‘to 


K, =Bal(hat+h.)*(1+12)/2, (81) 

from which it follows at once that 
ay = By = Bal(hath.)*/2, (82) 
V9 =Ca/Vq =Cal(he+ha)/2 J+. (83) 


For low attenuation h,., haX<1, consequently 
the phase velocity vs, exceeds the free space 
velocity c, which in turn is greater than the 
velocity cq in the dielectric d. For perfect media 
the phase velocity becomes infinite, and the 
attenuation and phase shift both vanish. 


Case III. Attenuation, v?>1: Equation (71) 
now becomes 


' K,* =Ba*(v?—1)(1+2h’), (84) 
with 
h' =(h.v? +ha)/(v?—1). 
Similarly it follows at once that 
a =Ba(v?—1)4f(h’); By=Ba(v?—1)!g(h’), (85) 


and for good conductors and dielectrics h’?/8<1, 
so that 


a, =Ba(v?—1)'; By=Ba(v?—1)*h’/2. (86) 


Here the attenuation is large and the phase shift 
small; in fact, for perfect conductors the wave- 
length \, in the pipe is infinite, and there is no 
wave phenomenon at all. The cut-off is sharp, 
unlike the gradual cut-off for the case of im- 
perfect media. 

For good conductors a,- may be determined® 
by integrating the normal component of the com- 
plex Poynting vector’ S = (EX H*) /2 over a closed 


18 Reference 7, pp. 135-137. 
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surface whose ends are planes perpendicular to 
the axis of the pipe and whose sides are the inner 
pipe walls. The approximation made requires that 
the surface current of the good conductor be con- 
sidered essentially the same as that for a perfect 
one, and is in accord with the assumption that 
(h?/8)<«1. Many different (but equivalent) ex- 
pressions for a,- have been obtained,®"’ and so 
we will content ourselves with listing only two 
forms here, to which we shall have some recourse 
in Part Il for the purposes of calculation. A 
detailed discussion of the variation of a, with 
frequency and the cross-sectional dimensions of 
the tube is beyond the scope of this paper. 

For the H modes of a rectangular pipe of 
width a and height b we have 


(=) a 2) } 
es ——3 
Qg b me 


fee ee, 


where p=a/b and n, and ny, are the number of 
half-wave-lengths in the x and y directions. The 
quantities €,, u-, and ¢, are the absolute permit- 
tivity, permeability, and conductivity of the 
pipe walls, and r, is the surface impedance 
obtained from the usual skin thickness formulae :” 


r,= (wu-/20-)}. (88) 


For the H modes in a circular pipe of radius a 
we have 


ot ke n? 
an= (=) (1 -»)-f 2 +e (89) 
@ \ be (AY, (n, »))?—2* 


where, following the notation of (52), m is the 
order of the Bessel function and p is the number 
of the root of the equation from which (ynp)reat 
is obtained. The attenuation aga is found at once 
from (68) and (78). 





THE CHARACTERISTIC IMPEDANCES 
V 


Unlike the conventional line, there are two 
possible definitions of characteristic impedance 
for the wave guide, one for the E and H modes, 
respectively. Both expressions are obtained by 


2° Reference 7, Chapter IX, Sec. 9: 17 (46). 
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formal analogy with the ordinary theory. Before 
discussing the character of this analogy let us 
examine the analytical behavior of these im- 
pedances. From (23) and (24) we have 





} me = K,“?/ Y.=R..(1 —1 ce); (90) 
Lem = LZm/ Kg = Rem(1—tdem)- (91) 
Dropping the subscript g we have from the 
previous section “ 
BY +qOfg@ 
(q=1) Ree= (1+ hg?) * 
WE 7 
: (92) 
a) — BOA 
BD gD fg” 
and go 
wuss 
~ stasiad ’ = a (2) /Q(2) 
(q=2) Rem o(2)24 g(2)2” Gem =a /B, — (93) 
and 
Xe{h} = —de{e}Re{h}, (94) 


where the upper line applies as g=1 and the 
lower for g=2. 

The H modes yield simpler expressions for the 
characteristic impedance, since the dependence 
on the imperfect dielectric 4a does not appear 
explicitly in Rem Or @em and occurs only in ay. 
This follows from the fact that for the H modes 
there is no primary axial E field (although a 
small one does exist in the instance of imperfect 
conductors and dielectrics), and consequently 
there is no appreciable displacement of free 
charge in the direction of the tube’s axis, a 
simpler state of affairs than for the E modes, 
which have both longitudinal and transverse 
dielectric currents. From the dependency of 
Z-{%} on K,™ we must then distinguish three 
cases. 


Case I. Transmission, v?<1: The general 
relation for Z.. and Z.m may be obtained from 
(73) and (90), (91). We shall consider here only 
the case of good conductors and dielectrics, i.e., 
eK, ha«K1. We have 


i 
R..= (*) (1—v2)!s bco= al /BO —ha™, 


€ 


be<Ree, (95) 
and 


Mu 


4 
Rem = (*) (1-— v2); bem =a /B™, 


€ 


bn<Rim. (96) 
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For perfect conductors and dielectrics ¢,.. and 
dem Vanish. 
Case II. Cut-Off, v?=1: From (81) we have 


b 
Reem | “(he +h 2| (1+ha)/(1+ha2) 
€ 


7 
=| +h) a], (97) 


€ 
Gee =(1— ha) /(1+ha), (98) 
be -} 
| -| ~(ha® +h.) 2| - bem = 1. (99) 
€ 7 


For perfect conductors and dielectrics R.» be- 
comes infinite, R.. vanishes, and ¢.. equals 


dem=1. 


Case III. Attenuation, v?>1: Here we have 
from (84)—(86) 


Tp fA’ 
R.e =| —(v? — 1) —+h, }; 
LE 2 


PA)’ 1 
; +h? ’ 
- 


(100) 
Dee = 





and 


m } 62’ ad 
Rem = rai ” (y?— 1) “23 Qem = (h® / 2) 7 (101) 
€ 2 

and for perfect materials R..=R.m=0, and 
duom=Oce= ©. The reactances X,.}%,} may be 
obtained from (92)—(101) with the aid of (94). 
The above results are summarized briefly in 
Table I.7! 





TABLE I. 
¥ Zee Zem 
Trans: <1 =R.. = Rem 
Cut-off: »? = = +0 
Atten: vr? >1 4X ..<0 4X. >0 
CONCLUSIONS 
VI 


As we have seen in the previous work, there 
is a formal analogy between potential and stream 
functions and the impedances of the wave guide 
theory with the potential difference V, total 


1S. A. Schelkunoff, Proc. I. R. E. 25, 1457 (1937). 
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current J, and impedances as given in ordi- 
nary transmission line work. The wave guide 
problem is three-dimensional, whereas the con- 
ventional line is essentially one-dimensional. 
Since it is always possible to separate the axial 
variable z from the transverse variables u, v, the 
correspondence between the conventional theory 
and that discussed here for hollow pipes of 
various cross sections is seen to lie in the fact 
that the dependence on the axial variable z is 
the same in both cases, and thus the properties 
of one go over without change into the properties 
of the other provided they are functions of z 
only. The propagation constant K“, however, 
depends in part on the transverse nature of the 
transmission circuit and on the mode excited, so 
that in general, unless we specify a particular 
mode and a particular cross section of pipe, the 
analogy between the two theories is only formal. 

This formal equivalence is not physically 
complete, for as we see from (18), (36), (37), 
(51), and (52), the stream functions do not 
reduce to the total current J in the limit of very 
small cross sections since 7(u,v) does not then 
reduce to unity. The total currents are still 
defined in the dielectric between the conducting 
walls, and the surface currents must be obtained 
from them by differentiation. It may easily be 
shown from (36) or (37) that when the circuit is 
infinitely long, the quantities Z.|%.} are correctly 
defined (just as V/J for an infinite line yields Z, 
for the line), for a given mode and a given pipe. 
In this case we may reasonably speak of an 
attenuation constant a,, phase shift 8,, wave- 
length in the guide \,, and ®{® and p®. Their 
physical interpretation is the same as for the 
corresponding quantities of the conventional line. 
Thus the familiar /, c, g, r (loop inductance, 
capacitance, conductance, and resistance) of the 
line may have their counterparts in the instance 
of the hollow pipe, but subject to the cross- 
sectional dependence mentioned above. When 
the termination is a device such that p,= «, the 
circuit, characterized by 


as before, will be 


traveling waves and a standing wave ratio of 
unity. The characteristic impedances Z,{%,} for 
all coaxial or hollow tubes reduce to an expression 
which is independent of the transverse dimen- 
sions of the pipe, and is (u/e)!, which has the 
value 376.7 ohms if the dielectric is air. This 
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means that for matching purposes it is not 
sufficient to equate Z.{%,} of the two guides (as 
it was for two lines), since for tubes of different 
cross sections and excited in different modes the 
driven tube will quite obviously not have a unit 
standing wave ratio. Matching may be achieved, 
as before, when the standing wave ratio is one, 
but to accomplish this will require adjustments 
of the line variables which are not contained in 
Z.{%} alone. 

In a similar fashion the terminal impedances 
Z)”, Z,® are physically analogous to the 
terminal impedances defined for the conventional 
theory,!:"! but only for a particular mode and 
pipe. The definitions of impedance (28), (30), 
(33), and (35) are not, however, generalizations 
of Zo, Z, of the ordinary two-wire or coaxial line, 
as explained above, since Zo and Z, do not 
reduce in the limiting case to Z» or Z,, but are 
dependent on the distribution of current within 
the pipe, i.e., on the stream functions, and not 
on the total current, as are Zp and Z,. The 
total current in these cases is of little significance 
in characterizing the behavior of the guide. 

We have already developed in some detail the 
formal reduction of the wave guide problem to 
that of the conventional line. From the stand- 


point of practice the most important results are 
that— 

(1) the formulae of the ordinary theory are 
applicable to wave guides and cavity resonators, 
subject to the restrictions and interpretations 
discussed above ; 

(2) since it is possible to determine p,“ and 
®,™ experimentally from any termination (as 
will be mentioned briefly in Part 11), the general 
relations of the usual theory may then be used 
to obtain specific experimental and theoretical 
answers to an almost unrestricted number of 
transmission circuit problems involving wave 
guides and resonators. Our success in obtaining 
analytical expressions for various terminations 
will depend on our ability to determine the field 
distributions, which is an extremely difficult 
problem for even the simplest geometries, and 
can be successfully treated in but a few cases. 
However, meaningful results may always be 
obtained experimentally in terms of p,“” and ®,’. 

In Part II we shall consider a few practical 
applications of the theory, with several experi- 
mental examples and an analytical evaluation 
of the terminal and input impedances of a closed 
guide or resonator containing several sections 
of imperfect dielectric material at one end. 
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In Part II a number of practical applications of the theory outlined in Part I are considered. 
Formulae for power transfer, efficiency of transmission, standing wave ratios, and insulating 
disks are given, and the experimental determination of the line and terminal parameters is 
treated briefly. A number of examples illustrating the theory are discussed: (1) a rectangular 
wave guide terminated by a section of pipe containing a horizontal slit at one end and a flat 
horn of 45° aperture at the other, (2) the same guide terminated by a section of tube with 
a horizontal slit at one end and opening into space at the other, (3) the rectangular guide 
with a variable horizontal and vertical slit, in either case followed by a closed quarter-wave 
section of pipe, (4) a closed cylindrical tube containing several transverse slabs of imperfect 
dielectric material. In (1)—(3) the terminal functions p and ® are determined experimentally, 
and in (4) analytical expressions for the terminal and input impedance are obtained. Experi- 
mental and theoretical curves illustrating the examples are included. 


VII 


PACE does not permit a detailed discussion 
of the many applications of the formulae of 
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Part I to the practical problems involving trans- 
mission circuits. The actual derivations used in 
the applications outlined below follow immedi- 
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ately from the theory of the previous sections.* 
The development given in the above reference 
pertains primarily only to parallel and coaxial 
lines of transverse dimensions small compared 
with a wave-length, but, as indicated in Part I, 
these results may also be applied, after some 
modification, 
guides and 


to circuits consisting of wave 
cavity resonators. Some of the 
formulae are listed here for convenience. 


(a) The Transfer of Power and Efficiency 
of Transmission 


A wave guide of length s has a load impedance 
Z, given by (40). The power supplied at the 
input end of the circuit is 


Ay” 2 


Wo | 2] 


Py =| q=1, 2, (102) 
where A," and ¥, are given by (37) and (36) 
suitably specialized at z=0. The power supplied 
to the load is 

, g=1,2, (103) 
with the help again of (36) and (37), at z=s. 
We have finally 


(=) sinh 2p,” —@,}%,} sin 26,‘ 
Py sinh 2(a‘?s+ p,”) 


— de} in} sin 2(B°Ys+4,) 


(104) 


From this it is a simple matter to compute 
(P,/P >). In all practical cases of large loads 
p.™ is sufficiently large so that we may neglect 


terms in ¢-,}%,}, and (104) becomes 


sinh 2p,‘” ” 
(P,0/P0)) mS 20 __ (gs) 
sinh 2(a‘s+p,‘”) 


provided 


sinh 2p,‘ >@-}%,}. (106) 
For maximum efficiency we find 
p\V= a, 
and 
P,&? Py? =exp ( —2a‘s), (107) 





* The reader is also referred to R. King, J. App. Phys. 
14, 593-599 (1943), for a more complete consideration of 
such topics as power transfer, efficiency of transmission, 
standing wave ratios, impedance matching, insulating 
stubs, and the experimental determination of the line and 
terminal parameters. 
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(b) Standing Wave Ratios 
or ° ° ~(a) (q) 
rhe standing wave ratios Stream, Spotential, May 
be defined as follows, by analogy with. the con- 
ventional line theory, in terms of the stream and 
potential function ratios: 
SY = Ay | max/ | A 1“) min; 


ae iS gigghciiot (108) 
Se —— 7 V; “| max Wy “| mings 
and 
Sp‘ =| vy, | max/ | S Pe | min; 
Sp‘ = A hi max / A 1?) | min- 
These ratios are equivalent to a comparison of 
|\£ | max E mings H max iT min OF any quantities 
whose maximum and minimum values may be 


(109) 


measured by a properly inserted loop or straight- 
wire probe, which is sufficiently loosely coupled 
to the circuit so that it does not disturb the dis- 
tribution of the quantities measured. Unless 
otherwise specified, it will be understood that 
the distance w=s—z is to be kept as small as 


‘ () . = 
possible, to keep S;,, single valued. Using (60) 
and (61) maximized and minimized with respect 
to w we obtain 


(q) 
cosh (a‘? Wimax + ps” 


) 
So = = So”, (110) 


a (a)en)°?? (q) 
sinh (a? wii, + ps'? 


where for S,“ and S, with n integral 


(q) 


Wmin =MAg/2—8,° /B™, (111) 
2n+1 
Wen = ( ae ee Bo”, (112) 


with 2 chosen (independently) as small as pos- 
sible in each case. For S,‘” we must interchange 
the subscripts max. and min. in (111), (112), and 
similarly for S,“. For low loss tubes, i.e., pipes 
with low over-all attenuation which are used to 


and 


. ° (q) (q) 
drive a power-absorbing load, Wmin and Wmax are 
always small enough so that we can write 


a) 
(a wirin, max) <ps, (113) 


and 
(a) 
S,'p=coth ps”, (114) 


a . ‘ (q). 
which determines p,‘” if S,,, is known. 


(c) Impedance Matching and Insulating Disks 


The problem of matching two sections of a 
guide, or inserting a section of tube between a 
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load impedance Z, and a long pipe cannot be 
solved in general by any of the analysis of Part 
I, except in the elementary case where the two 
guides have identical cross sections and are 
excited in the same mode. As was pointed out 
in Section VI, the characteristic impedances 
Z.\%} are not capable of the analogous physical 
interpretation of conventional transmission line 
theory when the pipes are dissimilar, and conse- 
quently one cannot equate, for example, Z.{%! 
of one section to the input impedance of a given 
length of pipe with a load Z,™ and obtain 
satisfactory results. Matching may be achieved 
experimentally, however. Consider again, for 
example, the problem of matching a rectangular 
guide supplying power to a circular guide excited 
in the same mode and terminated in an arbitrary 
fashion. One simply determines A, and F, 
for the terminated pipe; then we may treat this 
entire pipe as the termination of the other, 
whereupon A, and F, for the circular guide 
become p,“ and ®, for the rectangular guide. 
With suitable slits in the rectangular tube to 
vary the reactance and a sliding section of pipe 
to vary the resistance (as well as the reactance), 
it is possible to obtain a match, which as before 
is characterized by a standing wave ratio of 
unity. 

At high frequencies it is often very convenient 
to be able to use high impedance sections of 
tubes. Such sections are easily obtained, through 
the use of a metal piston at the end of the guide. 
If the over-all attenuation [or A,? and ¢.{%}?] 
is small compared to unity, as is generally the 
case in practical circuits, the maximum input 
impedance of a section of length s is 


(Rin) max = Retin} /La%st+p.], (115) 
which occurs at 
B%s=nr—O,™, (116) 


Thus the highest resistance for a given guide 
occurs when the over-all attenuation of the 
section is smallest, i.e., when the denominator 
of (115) is reduced as much as possible. The 


conditions for a piston are then from (43) and 
(44) 


p\%=0; 6,%=7/2, (117) 
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where, for »=1, we have 
(Rin) max = Refi} /as’, (118) 
s’ = 4/28 =), /4, (119) 


where 


with A,“ the wave-length in the pipe. 


(d) Experimental Determination of the Line and 
Termination Parameters 


For a more detailed discussion of this topic 
see Part I, reference 1.* We list here only the 
relevant formulae, with a brief description : 


Aj =2(se—s1), (120) 


where s; and sz represent the location of two 
successive field maxima [or maxima of A; or 
Vv, ], provided the width of the resonance peaks 
are small so that (65) is true. For the measure- 
ment of the attenuation constant of the circuit 


we have 
a® =B (As, —As,)/2(S,—51), (121) 


where As, and As, are the widths of two resonance 
peaks at 0.707 of maximum and s,—s; is the 
distance between them. The position of the de- 
tector and source must not be changed during 
successive measurements, nor must the degree of 
coupling be altered. 


’ , (q) ° ° 
The phase functions ®,,9 of a termination 


are determined readily from (62) by 
(q) 


, 9 =7/2+8M(si—Sia), (122) 
where s; is the resonant length of the guide 
terminated by a piston and sie the resonant 
length when the unknown termination is added. 
Condition (63) must be satisfied in this case. 

The attenuation function is obtained from 


ps =B (Asy,—As1)/2+a™(si—Sip), (123) 


where As; is the width of the first resonance 
curve with the metal piston terminating the 


guide (p.e%0), and As, is the width of the first 
resonance curve with the termination of unknown 
p:; s; and sy, are the respective resonant 
lengths of the two cases. As before it is very 
important that the degree of coupling be un- 
changed during the measurement, as well as the 
location of the source and detector with respect 
to the end at z=0. 


* Part I, reference 1, pp. 598-599. 
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Fic. 1. Resonance curves obtained for a rectangular 


brass pipe excited in the H;,; mode and terminated in a 
horizontal slit 1.2 cm wide, with 11.5 cm of pipe extending 
beyond to an open end. The termination is capacitative 
and is characterized by p, =0.54 and ®,® =0.79 rad. 


The Q of a resonant guide or cavity may be 
defined in the customary way by 
Qo =f, /26fM =s,/As,, (124) 
provided (65) is satisfied or (6f)*<f,?. Here df is 
the change in frequency required to reduce the 


resonant amplitude of the fields to 2—! of their 
maximum at f,. From (64) we have 


Ow = Bs, 2(a 5,9 + po +p,”). (125) 
. (q) 
For po.» =0, we have 
Q@ =B /2a, (126) 


(e) Examples 


Let us consider a number of practical examples 
which are not easily handled by analytical 
methods but which are immediately capable of 
experimental treatment. 

To illustrate the case when p,™ is small we 
have a long rectangular guide excited in the 1, 
mode and containing a source and detector at 
one end and a horizontal slot 1.2 cm wide at the 
other, followed by 11.5 cm of pipe leading to an 
open end (see Fig. 1). The guide contains a 
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calibrated moveable piston at the end z=0. The 
source and detector are placed so as always to 
remain at the same distance from this end and 
from each other. This is important, as otherwise 
the attenuation and phase shift due to the source 
and the detector, which are contained in the 
terminal functions po and 9, would vary 
with successive measurements and would not 
subtract out, as they do in (123). We would then 
be unable to determine p,“ and ®,™ accurately. 

With pistons at both ends the attenuation of 
the circuit was measured with the help of (121) 
and found to be negligible, as was p,® for the 
piston. From (120) and (123) it was found for 
the rectangular guide of width a=6.35 cm and 
height 6= 3.15 cm that 


A, =90.150m; Asn; =0.006 m; 


B, =22r/d, = 41.9 rad/m; (127) 


po’? =B,Asn1/2 =41.9(0.006) /2 =0.125. 


The termination to be tested was added and 
another resonance curve obtained from which it 
was found that 


Asne=0.0315 m and sny—Sn2= —0.0185 m, (128) 
and from (122) and (123) we have 


ps‘? = 41.9(0.0185) /2 =0.54; 
(129) 
®,'? = w/2—41.9(0.0185) /2 =0.79 radian. 


Since ®,° <2/2 the termination is capacitative, 
as can be seen from Fig. 6, reference 1, or (42). 
Figure 1 is a plot of the data from which these 
terminal functions were obtained. 

When p,“ is particularly large, 1 or greater, 
we cannot use (123), as (65) is not satisfied. 
For practical purposes the guide is then nearly 
equivalent to one terminated in its characteristic 
impedance. Values of ®,“ are not in this case 
required. To determine p,“ we measure the 
standing wave ratio and use (114). For the 
guide of the previous example, excited in the Hu 
mode and terminated by a horizontal slit fol- 
lowed by a telescoping section of line which opens 
into a flat horn of 45° aperture, the standing 
wave ratio was found to be 1.2, and from (114) 
we have 

ps =coth-! (1.2) =1.20. (130) 
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A third example consists of determining , 
for the same rectangular guide terminated first 
in an adjustable horizontal slit followed by a 
quarter-wave-length section of closed pipe, and 
then in a variable vertical slit, also followed by 
a closed quarter-wave section. The measure- 
ments were made in a fashion similar to the first 
example, with the help of (122), and in each 
case ps“ and as were found to be negligible. 
The phase shift @,@ is shown in Fig. 2 as a 
function of slit width. 


CYLINDRICAL WAVE GUIDES AND CAVITY 
RESONATORS CONTAINING SEVERAL 
IMPERFECT DIELECTRIC MEDIA 


We will now determine the terminal and input 
impedance of a cylindrical wave guide containing 
several sections of imperfect dielectric materials 
at one end of the closed pipe, for the general case 
of forced oscillations from a _ closely-coupled 
source. The analysis is equally applicable to the 
case of loose coupling, where we may consider 
the region containing the source to be of negli- 

















Ls 
of Eee ee 
zane ® 
(2) 
$; 
) 
open s closed 
(2) Horizontal slit 
T 
$ 
TD, 
(2) 
4 
0 
open s closed 


(b) Vertical slit 


Fic. 2. The terminal functions #,® for a rectangular 
guide terminated in (a), a horizontal slit followed by a 
closed quarter-wave section of pipe, and (6), a vertical 
slit, followed by a closed quarter-wave section. 
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Fic. 3. Schematic representation of a wave guide or 
resonator containing a source and detector at one end and 
terminated at the other by several sections of imperfect 
dielectric material. 


gible longitudinal extent (cf. Section III), i-e., a 
plane distribution, and then we may quite 
reasonably treat the pipe as a cavity resonator, 
when, of course, it is excited in one of its natural 
modes, either by an adjustment of the termina- 
tion, the exciting frequency, or the length of the 
tube. In the cases of loose coupling the 
resonance condition is given by (62). 

The circuit with its various dielectrics is shown 
schematically in Fig. 3, where the generator 
and detector are confined to a fairly small 
region at one end of the pipe, and the dielectrics 
are placed at the other. Another practical 
arrangement is to have the source and detector 
somewhere about the middle of the tube, and a 
slight modification of the analysis (see Section 
III) enables us to treat this case at once. 

Although we are dealing with forced oscilla- 
tions, the wave equations (11) and (12), speci- 
fying the fields in media one-three, are homo- 
geneous, since the regions under consideration 
contain no sources (or sinks). The propagation 
constants for the various media (j=1, 2, 3) are 
given by Eq. (20). (See also Section IV.) We 
have remarked before (Sections III, IV) that the 
propagation constant K“ depends on the geom- 
etry of the tube and on the type of mode excited, 
as well as on the electrical constants of the 
medium (k,”) and the frequency of the source. 
In the following analysis we will assume that 
the pipe walls are perfectly conducting, and 
only the dielectrics are imperfect. Since the 
contribution to the attenuation a,4 due to the 
dielectric may be added directly to that due to 
the metal walls, viz. (78), for imperfect but 
highly conducting metal, the rigourousness of our 
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theory is not materially affected in the practical 
case.* The development of Part | is quite general 
for all cross sections for which (18) is soluble. 
Here we will consider specifically only circular 


pipes. 

The fields are expressible in terms of the 
vector stream and potential functions AM by 
virtue of (5), (8), (10), (14)-(16). For the E 
modes (¢=1) we have 

FA 184 | : 
E) =| i,;— +i. — ig (D2 /¥. 
Ordz r d0dz 
(131) 
OA || ; 0A," - ; 
H‘) =i, —1|y ; YF.=twete. (132) 
roe or 


Note that because of our assumption of perfect 
is real; 


walls, y‘” see (57). For the HT modes 


109 





(2 Ra) ——t —— ~— 
S,(in)max $,Un)max-——— 
J 


_— . 
+X. (inymax -X5 (in)max 
—V¥V,= 58 —t— 

Vs = 803 —— 
— fg= 24£,=24€,- 
AK =10-0cm 
= a = 365-10 ‘m 
=, = 240%m 
Hy 


Terminal impedance at X, 


Fic. 4 = 16.8 cm for a section 
of cylindric al brass pipe of radius a=3.65 cm, brass 
= 1,43 107 mhos/m, excited in the H,; mode and containing 


a thickness /2 of dielectric (¢.=2-10~* mhos/m, ¢€=2.4¢), 
as a function of J3. Here a; =a;=2-11-1073, a2®@ 
= 2.93-10°3 nepers/m, 82°? =83.4 rad/m, 8: =83° =37.5 


rad/m. This also represents the input impedance of a 
— or resonator Containing a section of dielectric at one 
end (/;=0) as a function of the tube length s =/5. 


* Note: To date there has been no completely rigorous 
treatment of the finite, imperfectly conducting cylinder. 
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.medium one, (j= 1), 





(q=2) the fields become 








0A aA, | 
E(?) = —j,—_ + i,—_, (133) 
r0e0 or 
0°A, (2) A. (2) 
* ards  rabaz 
; , (134) 
Zun=twy, (t= —1) 


The quantities i,, i2, is are unit vectors in the 
r, 6, z directions, respectively (see Fig. 3). We 
have then to satisfy the boundary conditions 
(53), (54) at the tube walls (r=a), and at the 
dielectric interfaces, 2=1/), 1; +/2, 1; +1l2+/1s, where 
it can be shown with the help of the equation 
of continuity that the conditions (54) become 


(q) ° 
n-¢;*E;?=n- ees 1Ejs1, €*=€—10/w, 


(135) 


ee 


nXH;? =nXH;, 7=1,2. 


~ 


The surface charge and current densities on the 


cylindrical walls r=a are given at once by (54). 
We find for both modes that 
T(u,v)=Jn(ynpr) cosn@, g=1,2, (136) 


is a sufficient solution for (18), where p indicates 
(q) 


the pth root of the equations determining y,,, 
i.€., 
(1) 


Tn(y¥np@)=0; Jn'(ynpa)=0. (137) 


Typical solutions of (11) then take the form for 
(OS25/),=s): 


Ar, 0 5)= (C0 3 | 
lsinh K,z | 


| yr (r, 0), (138) 


| cosh 
A: 


V(r, 0, 2)=( gV 


\sinh 


and for medium three, j7=3, (J) +/.S25)4+/2 
+1];=50): 

| cosh K3°?(so—sz)} 
A;(r, 6, 2)=C;“ e | Tq V(r, 6). 

| sinh K;° (so—z :) | (140) 
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Fic. 5. The p,® function for the termination of Fig. 4. 


This also represents A,® for the hollow pipe of length /; 


terminated in a dielectric section /». 


The complete solutions, of course, contain the 


time factor e* in product with (138)—(140). 
Here the upper and lower lines apply as g=1, 2, 
respectively. From the boundary conditions (53) 
and (135) at the various dielectric interfaces, 
and solving for Ci/By and Di®/Byo@ with 
the help of (30), (34), (138)-(140), we obtain 
finally for medium one, g=1: 


A,(r, 0,2) =T(r, 0) 
jr sinh Kyo (s—2) +cosh Ky) (s—2 z) | 


——_____—_— +, (141) 
| rey sinh K,s-+cosh Ks | 





and for the H modes we have 
A,(r, 0, 2) =T®(r, 6) 
{sinh K,(s—2z )+r cosh K,‘?)(s s—z)| 


- ; —+, (142) 
sinh K,@ Yst+T® cosh K,s 





where 


tanh K,“/,.+7. tanh K3%ls 
‘(q —_—, (q) 
Tq) = 71 





1+72 tanh Kl, tanh Kl; 


(143) 
=f /7,@ (143a) 
and , . 

7°) = €o* Ky) e:\*Ko"); 

: (144) 
ro) = @o*K;()) e;*K.”, 
71 = wi Ko /pok; 

(145) 


72 = usK wok; 


By applying the appropriate analysis (26)—(35), 
we may obtain (36) and (37) and expressions for 
Z, in terms of T®, where T™ is given by 
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(136) and A,.™, VW. by (27) and (34). Media 
two and three form our termination, as can 
readily be seen from Fig. 3, and their field 
distributions may be found from (139) and (140). 
The relations between the various amplitudes 
Ci, C2, C3, Dy, etc., may also be deter- 
mined from the boundary conditions (see Fig. 4). 


TERMINAL AND INPUT IMPEDANCE 


The terminal impedances of a wave guide 
containing several sections of dielectric material 
at one end become at once, from (23), (24), (31), 
(35), (143), and (143a) applied to (141) and 
(142), 

—iK /wes*) — 

z=} : _ lr, g=1,2, (146) 
icope/ Ko 

where the upper expression in the brackets is to 

be used when g=1 and the lower when g=2. 

From the general relations (39) and (40), p,‘® 

and ®, of our termination are defined and may 


be obtained from (43) and (44). In a similar 


fashion zs may be obtained from (49) and (50). 
(See Figs. 5 and 6.) 
When the dielectrics are perfect, we have 





K{? =1 fo, 7=1, 2,3 ’ : 
j B; (7 ) (147) 
72°) = €083") / e382; — 9? = 382" /poB3. 
Then (146) becomes 
1B) / 
Z,” = g wea| 
iws/ Bo? | 
tan B2%l.+72 tan B;l 
| 2\V bo 3Vl3 . (148) 
1— TT? tan BoM, tan B3Mls 
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Fic. 6. The #,® function for the termination of Fig. 4. 
This also represents F,® for the hollow pipe of length /; 
terminated in a dielectric section /2. 
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Equation (148) is a pure reactance, as one 
would expect, since there are now no dissipative 
elements in the termination when the dielectric 
conductivities are zero. However, the general 
expression (146) is complex and has a real or 
resistive term R,“ which is positive and non- 
vanishing. When 


A,M€1, o-}%}2<1, (149) 


which is characteristic of a guide with low over- 
all attenuation, low attenuation per unit length, 
and negligible leakage conductance, the extreme 
values of R, and X, become, upon simplifi- 
cation of (41), (42), (49), (50), 


(R,™) min = Rel ajAs. 2 (150) 

at 
tan (Bs %ls+,"») =6;A,'s a;”, (150a) 
(Remax = Rel in} /As.2 (151) 

at 
8;%l;=n4—9, », (15la) 
(X,) max = LRe{%} / 2A,'o](1 — a3? B;), (152) 


at 
) 
B3%l3 = (ne —$,'2)(1—as, BS?) — ont (152a) 


10* 


To 0 
| 
}|| 
mhos/meter 
vo 4 1 





2 46 


Fic. 7. The terminal impedance at 4,=16.8 cm, Hi; 
mode, of a section of the pipe of Fig. 4 containing a 2-cm 
length of dielectric, as a function of dielectric conductivity 
C2, l; == (), 
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Fic. 8. p,® for the termination of Fig. 7 as a function of 
dielectric conductivity oe. 


> (q) 
(X PS ) min = [R. ; i ’ 2A 8, 2 | 


X (1 +a3°7/83), (153) 
at 
Bs = (nx—,3)( 1+a3”?/B3;) +pns (153a) 
where 
(q) (q) - 
A ,.2=a3M13+ ps, 2. (154) 


The functions ~s and ©," are obtained from 
(43), (44), and (146), with /;=0. For the guide 
terminated in Z,, (146), we may simply obtain 
the extreme values of input resistance and 
reactance by replacing the subscript s on R,‘® 
and X, by in and changing /; to s in the above 
formulae, as well as a3, 83; to a;, By, 


respectively. The quantity F ac4 becomes A,, 
and is given by (50), where p,“ and #,@ are 
obtained from (43) and (44), for the impedance 
of the two dielectric termination (146). Figure 4 
illustrates the variation of this impedance with /3, 
when medium three is air. As can be seen from 
Fig. 4, and (151)—(153), the terminal (or input) 
resistances and reactances approach critical 
values of approximately half a megohm, which 
are the same order of magnitude as the values 
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Fic. 9. &,® for the termination of Fig. 7 as a function of 
the dielectric conductivity oe. 


obtained in the case of the conventional two-wire 
or coaxial line.* Figures 5 and 6 are calculated 
from (43) and (44) and give the variation of the 
terminal functions p,® and ®,@ with 13. The 
phase function ®, varies between r—6 and 6 
(6=10-*) and changes very abruptly in the 
immediate neighborhood of anti-resonance. This 
change would be discontinuous for perfect media, 
and p,;® would vanish. For the low loss termi- 
nation or circuit, p,® (Fig. 5) varies over a 
fairly wide range and reaches a maximum value 
at anti-resonance, but is never at any time 
greater than about 1.5-10-* in this instance (see 
also Fig. 7). 

For the important special case of a one- 
medium termination (146) reduces to 


—1iK, wer™| 


Zz, =| | tanh Kos, l1,;=0. (155) 


icoue/ Ko 
When the dielectrics are perfect this becomes 


1B2'” /wes 


Z,” = 





; tan Bos, (156) 
twp /B\ 

which like (148) is a pure reactance. Figure 7 
illustrates the variation of Z,® of the one 
dielectric case (155) as the conductivity of the 
medium is varied.' Antiresonance effects are 
noticeable when ¢2 reaches a certain value; a 
brief inspection of (155) indicates that this is 
quite possible, for instead of changing /; (or /2), 


*See Part I, reference 1, Figs. 2 and 3. 
1 Kennelley’s Chart Atlas is helpful when dealing with 
hyperbolic functions of complex variables. 
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as was done for Fig. 4, K2™ is now varied. 
Figures 8 and 9 show p,” and ®,” as functions 
of o2. When o2>0.18 mho/m, e=2.4€, the 
termination is inductive and in effect lengthens 
the guide, though not greatly; for o.<0.18 it is 
capacitative and shortens the pipe. Figure 10 
illustrates the effect of conductivity on the 
attenuation and phase shift constants of the 
medium [a constant attenuation for the metal 
guide is included, additively when (77) permits; 
otherwise (76) is used]. The phase shift 6; is 
affected much less critically than a; when 
o2<10-! mho/m. 


BEHAVIOR OF THE GUIDE AS A RESONATOR? 
IX 


The theory of the previous sections applies in 
general for the hollow pipe as a guide or reso- 
nator. For resonator behavior the source is 
loosely coupled to the hollow cavity, and we may 
treat its distribution as being essentially plane, 


103 10? 
(2) 
Be (2) 
radians/ Qa, 


meter nepers/ 


meter 


1o* 
H,, mode 


Xr = 10cm 
f, = 2:10? m 
V =518 

€, = 24, 


+ 


6, mhos/meter 





1tt : ~ 
10 10 10 10 1" 1246 


10°? 
Fic. 10. The attenuation a2® and phase shift 62® of the 
pipe of Fig. 7 as a function of dielectric conductivity a2. 


2D. Middleton, Phys. Rev. 63, 343 (1943), gives a 
discussion of cylindrical cavity resonators containing 
several perfect dielectrics. 
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as was pointed out in Section III. If the coup- 
ling is sufficiently loose, the frequency of the 
source for which the hollow pipe becomes reso- 
nant, i.e., when standing waves are present in the 
interior, is the natural frequency of the cavity. 
Equation (62) yields the condition for resonance. 
The Q of the circuit,* consisting of source, de- 
tector, and closed hollow pipe, must be high 
enough so that (63) is satisfied. When the stand- 
ing wave patterns become very broad we are no 
longer able to determine the maximum ampli- 
tudes accurately, and in fact the term resonance 
becomes meaningless. Analytically &) and ®,™, 
as well as po and p,”, may be quite impossible 


* Part I, reference 10, Chapter 10. 





to determine, but they can always be obtained 
experimentally, as explained in Section VII (d) 
and illustrated in (e). The fields and charge and 
current distributions in the resonant case may 
then be determined accurately from (58) and 
(59), or their amplitudes from (60) and (61), 
with the help of (54), except in the immediate 
vicinity of the source and detector, as was 
pointed out in somewhat greater 


Section III. 


The writer wishes to express his appreciation 


detail in 


to Professor Ronold King for his encouragement 
and many valuable suggestions in the prepara- 
tion of this manuscript, and for the data of 
Figs. 1 and 2. 





An Approximate Representation of the Electromagnetic Field in the Vicinity of a 
Symmetrical Radiator 


CHARLES W. HARRISON, JR.* 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts 


(Received January 24, 1944) 


The actual quadrature components of current along a symmetrical center-driven antenna 
of non-vanishing radius are represented by simple distribution functions. The magnetic and 
electric vectors in the near zone are then determined by carrying out the calculus of the field 
theory involved. The resulting expressions, while more intricate than existing formulas based 
on the assumption that the current in any antenna is sinusoidally distributed, represent 
approximately the field conditions about antennas used in practice. 


HE electromagnetic field at points outside 
an antenna of half-length #4 and radius a, 
and oriented as shown in Fig. 1, may be con- 
veniently calculated from the vector potential 


given by 
g fr * c- 
A,= — [ I, dz’. (1) 
4n J-1 R 


Here the universal magnetic constant Il=47 
<10-7 henry per meter. J, is the complex 
current amplitude in amperes flowing in the 
element dz’. 8 is the propagation constant in 
radians per meter, and equals 27/X, where X is 
the wave-length. R is the distance from the 
point (r,@,z) where the vector potential is 
computed to the element dz’. That is, 
R=[(2-2')?+r°]}. (2) 
° Now at the School of Engineering, Princeton Uni- 


versity, Princeton, New Jersey, and the U.S. Naval Re- 
search Laboratory, Washington, D. C. 
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The electric and magnetic field vectors in 
empty space may be computed directly from (1) 
using the relations 


B=curl A (3) 
and 


E=1/jwIIA curl curl A= —j(c?/w) curl B. (4) 


In terms of the universal electric constant 
A=8.85 X10-" farad per meter, and the magnetic 
constant II, 


c=1/(IA)*=3X108 meters per second, (5) 
B=w/c. (6) 


The components of the curl in cylindrical 
coordinates involving the z component of A are: 


curl, A=(1/r)(yA./y9), 
curls A= —yA./yr, (7) 
curl, A=0. 
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Fic. 1. Illustrates orientation of antenna relative to the 
primed and unprimed coordinate systems. 


Since rotational symmetry obtains, yA ./7y@=0, 
and accordingly (3) becomes 


Be=curls A = —yA./yr. (8) 


The components of the curl in cylindrical 
coordinates involving the @ component of B are 


curl, B= —7yBo/v7z, 
curls B=0, (9) 
curl, B=(1/r)(y/yr)(rBa). 


Therefore, the radial component of the electric 
field vector is 
E,= (jc*/w)(yBo/7z) (10) 
and the tangential component is 
£2=(—je?/w)(1/r)(y/yr)(rBe). (11) 
From a theoretical point of view the calcula- 
tion of the field vectors in the vicinity of an 
antenna is a relatively simple problem, if the 
distribution of eurrent is known. However, the 
exceedingly intricate form of the actual current 
distribution in antennas of finite radius seems 
to preclude a direct determination of the field. 
An alternative approach, which is probably 
satisfactory in practice, is to represent the 
quadrature components of current along the 
antenna by simple distribution functions. This 
is actually possible to a good degree of approxi- 
mation. 
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For the component ofZcurrent in-phase with 
the driving potential Vo* one writes 


I (= Bz’ —cos ) 
° 1—cos Bh 3) 


and for the component of current in quadrature 
with Vo* the expression is 
sin B(h— |2’|) 
I,’ . . (13) 
sin Bh 

In (12) and (13), Io’ and Jo’ are the current 
components at the input terminals of the an- 
tenna. Their values may be ascertained directly 
from the relation 





Io= Vo Yoo, (14) 
where Jo is the complex current at the driving 
point, and Yoo is the input admittance of the 
antenna in question. 

Naturally there are intrinsic difficulties in any 
representation of the current except the correct 
one. The reader is referred to the literature for a 
discussion of this point.'? 

Substituting (12) in (1) 


IT )”’ 1 pt" e—i8(R-2’) 
Bag = \f dz’ 
4mr(1—cos Bh) | 20-1 R 


1 (th e—iB(R+2’) +h e—ipR 
+5f — a one ah R as} (15) 


—h h 








Using (8) one has 





—_— pe TT 9” 
4x(1—cos Gh) 
OY 278 (Reth+z) e~" 
x 3e+ 162— —du 
YY J ip(Ri-h+z) UU 


7B (R1i+h—z) er 
+ $e~ 1-— —dy 
Y’ JY i8(R2—-h—z) v 


+h Y eiBR 
—cos ph a as. (16) 


~ 





In (16) thar 

n 
Ri=[(2—h)?+r?}}, (17) 
Rr=[(st+h)2+r2}}. (18) 


Upon noting that 
7 w2 e—w e7w2 YWs enw yw 
see —dw =—- —-——-——_ (19) 
yr dui W We yr W1 «Yr 

and carrying out the required differentiation, 








‘Ronold King and Charles W. Harrison, Jr., Proc. 
I. R. E. 31, 548-567 (1943). 

* Charles W. Harrison,. Jr., and Ronold King, Proc. 
I. R. E. 31, 693-697 (1943). 
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Ilo" 


B,"’ ——— 
4x(1 — cos Bh)r 


e~ 7BR1 
x (jR: sin Bh+(h—z) cos Bh)— ; 
l 
e ibRe 
+(jR2 sin Bh+(h+z2) cos Bh)— es 


+h e~ iBR 
—r?® cos Bh (sf —dz’ 
-— RF 


wth e 7B8R | 
+ { a’) - (20 
a: R3 tf | ( ) 
The integral +h e~iBR 
T= { —dz’ (21) 
a R* 


with m=1 plays a major role in the solution of 
the problem of coupled antennas.* The evaluation 
of (21) may be readily carried out, using numer- 
ical methods. 


From (13), (1), and (8) 


jul’ 
———|2 cos Bhe-i8®o 
” dar sin Bh 


—eibR1 —e~ 1BR2} : 


(22) 
Ro=(2?+r*)!. (23) 

Equations (20) and (22), when appropriately 
combined in quadrature according to (14), com- 
pletely define the magnetic field Bg about a 
linear radiator of finite radius. 

The radial electric field is readily 
from (20) and (22) by applying (10): 

jlo” 


4ewAr(1 — cos s Bh) 
| ; r? cos Bh 
x (s(h+2) sin Bh+ R.? - 


B(h+2)? cos a 


i R, J R 


: r? cos Bh 
- (s0-2 sin Bh+ = . 


1" 


B(h—z)* aed -j8Ri 
—;- 4 = 


e~ iBR2 e BR, 
-—f COS Bh (> -5-) 
| R,' R;3 


e ibR1 e~ iBbR2 
af —__ — —_—__ J + @*( Cia 
+6"( Ri R: ) . 


- jSi@R,— CiaR.+- jsiaRs) |}; (24) 


obtained 


E fe 
“r 


~ 3 Ronold. King and Charles W. 
Phys. 15, 170 (1944). 


Harrison, Jr., J. App. 
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: jl’ (z—h) 
J oe ae ee —e~ 7BR1 
deel sin Bh\ R, 
Et) ss a OT sey . (25) 
Rz Ro 


Using (20) and (22) with (11), the tangential 
electric field is found to be: 
—_ i] aA 
=— = |(2 sin Bh 
4mrwA(1—cos Bh) 
(h—z)(1+ 78R;) cos a) 
Ri , Ri 


4? 














h+2)(1+ 7 R» } 
+(ssin PO tii od mene) 
R,? 
e 7bR2 ath — iBR 
x —cos i( 12 f — head 
R. Bh{ 728 ( 
+h e i8R 
+ (2+ ?r?) — 5, a2 
e/ —h . 


+h e~ 7IBR 
- (9 [- — mo 





+h e~ BR | 

+ - ~#)) ; (26 

. “a (26) 
jl’ | e~ 7BRo e~ 78R1 asi, | 

E,/ =——\ 2 cos Bh - = . (27) 
4ncA| Ro R, Ro 


The local radial field about an antenna of 
non-vanishing radius is given by (24) and (25) 
combined in quadrature in accordance with (14). 
The electric field in the z direction is given by a 
similar combination of (26) and (27). All field 
expressions apply. to the case of an isolated 
symmetrical center-driven antenna of half-length 
h and radius a, or to antennas of full length h 
and radius a erected vertically over a perfectly 
conducting half-space of infinite extent. 

The general complexity of the relations given 
here precludes a general graphical portrayal of 
field conditions about antennas of specified 
dimensions. 

The opinions contained herein are those of the 
writer, and do not reflect the views of the Navy 
Department or the Naval Service. 
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Theory of Membranes in Shape of Developable Surfaces 


LEON BESKIN 
Consolidated Vultee Aircraft Corporation, Allentown, Pennsylvania 


(Received May 1, 1944) 


HE theory of equilibrium developed here- 

after is an application to developable 
surfaces of the general’ equations of equilibrium 
of membranes; such theory disregards local 
bending in the thickness of the membrane. The 
following analysis deals with general develop- 
ables, exclusive of cones and cylinders. 


I. GEOMETRY OF DEVELOPABLE SURFACES 


A developable surface can be generated by the 
tangents to a skew curve, which is called the edge 
of regression of the developable. The parametric 
equations of the skew curve are taken under the 
form 


x1=filrr); X2 = fo(A1); x3 = f3(A1). (1) 


Thus, the equations of the developable can be 
taken under the form 


x1 = fi(di) +rofi'(A1), 
X2 = fo(A1) +A2fe’(A1), (2) 
X3 = f3(A1) +ofs’(A1); 


‘ indicates the derivation. When 


A, represents the length of the arc along the edge 
of regression, it is more convenient to replace \a, 
which is the length of the tangent between the 
point of contact with the edge of regression and 
the point under consideration, by \2.—),. Thus, 
A2 represents the distance between the point 
under consideration and a fixed point P on the 
edge of regression, distance taken successively 
along the tangent, and along the edge of regres- 
sion. Then: 


1 = fi(A1) +(A2— Aa) fir’ (Ax), 
X2 = fo(A1) + (A2—Ai) fo’ (As), (3) 
X3 = f3(A1) + (A2—A1) fs’ (Ai). 


With this representation, the curves \; are 
orthogonal to the generators 2, as the condition 
of orthogonality which is 


in which the 


OLn, Chu 


m Or 


OX2 
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becomes 
(Ao—A1) (fi' fir’ + fo’ fo!’ + fe’ fs’’) =0. 


This condition is easily verified since \, is the 
arc of the edge of regression, thus: 


fi? +fe? +f? =1, 


which, by derivation, gives the condition above. 
As the generators of a developable surface are 
lines of curvature, the lines \;, orthogonal to the 
generators, constitute the second set of lines of 
curvature. 
We call p the radius of curvature of the edge 
of regression 


1/p=(fi'?+ fo? +fs'")!. (4) 
The metric elements of the surface, i.e., the 
ratios: 


€,=d5s,:d\; €o =d5So:dd2 


in which s;, s2 are the arcs along the curves- 
coordinates, and also the direction cosines 
Li, L2, Ls of the normal to the surface are equal 
to: 


the Ae Al 
e=(fl? +f? +fs)(A2—A1) =——;; e = 1; 
p 
Li= — pl fo!’ fs’ — fs" fo’); 
Li= — pl fa fi! fifi’); (5) 
L3= —p(fi'' fo’ — fe" fi’). 
Also: “ 
Xm an? 7 
—— = (No—A1) fim” — fm 
Or" 
0 Xm Xm 
ei. phe 
OA10X2 Od? 


The geodesic curvatures of the curves Aj, A2 are, 
as it appears from Fig. 1, 


1 1 1 
fs —=0. (6) 


sees 
ri AeA) re 
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Fic. 1. 


The principal curvatures of the surface are: 


i _— f3'”"| 
fi’ fal" fal" | p? 
1 : Dh (fe! fa — fa'' f:') z fi’ fi’ fs! 


R’ (AoA) (Lf)! (X2—Ay) 
1 

It is to be noted that the curvature 1:R’ is 

simply related to the torsion T of the edge of 

regression 


. (8) 
R’ €) ry 


Il. EQUATIONS OF EQUILIBRIUM 


The general equations of equilibrium are 
written between the applied surface loads and 
the stress flows. The three components of the 
applied loads per unit surface are: Z, along the 
normal to the surface; X,, X2 along the tangents 
to the curves-coordinates. The three components 
of stress flow (i.e., stress multiplied by thickness) 
are m,, m2 normal flows along the curves-coordi- 
nates, #2 shear flow along the curves-coordinates. 
The three equations of equilibrium represent 
one, the equilibrium along the normal; the two 
others, the equilibrium in the tangent plane 
the two latter being merely the equations of 
equilibrium in polar coordinates, the pole being 
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taken on the edge of regression 


— + 
1 On, 1 O[ m12(A2—A1)? | : 
— —4+——__ ——___“-+x,=0, 
€: 9A; (A2—A;)? Ode 
(9) 
1 A[m2(A2—Ax1) | 1 On 2 
Ao At OX2 €1 Or, 
nN 
———-+X2=0. 
AeA 


It appears that m1. and m2 can be calculated by 
successive integration : 


ny= —ZR’, 


1 r2 
N12>= — — | (A2—A,)? 
(A2—A1)? Jas 


1 re 
-f | mead Xs 
AoA, Al 
Oni2 D(Ai) 
“— -Idd\o+ — 


Rae She 
pC’ (Ai) pC(A1) D(A) 


= 22° --—— — ' 
(A2—A1)?_ (A2—An)® Aa Ag 





Fic. 2. 
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The two successive integrations introduce two 
arbitrary functions C(A;) and D(\;) which are 
to be determined by the conditions at the edges 
of the surface. m12° and m.° are the values of 
ni. and m2 when the functions C and D are equal 
to zero, the limits of integration being those 
indicated. In the case when the applied loads 
are constant in value and direction, along each 
generator, or more generally in the case when 
the loads are functions of a power p of \2x—A,, 
i.e., of the distance to the edge of regression 
measured along the generator, the stress flows 
have the following values in which we omit C 
and D: 


Z=Z,(A2—A1)?;- X1=X1,(A2—A1)?; 
X2=Xop(A2—A1)?;— M1 =(A2—A1)?*'74; 
ZX")! 

Vi = 





SA” (fa! fal — fal fe’) 





pt+1 
N12 = (A2—A1)?* yo + (A2—A1)"pr1 ° 
pt+2 11) 
1 ” Ov; \ 
nz = ———j Xip + p— |; 
p+3 Or, 
N2 = (A2—Ay)?**veo+ (A2—A1)? de 
pt+l 
+ (A2—A1)?' p?»1—_; 
p+2 
= OV. p* Ov; 
ve= —Xop+r1—p—}; i2=pr2——— —. 
aay p+1 ar, 
If p=—1, all three stress flows are limited to 


their first factor, and if p=0, nz is limited to its 
two first factors. 

The most important cases correspond to p=0, 
in which case the surface loads are constant 
along each generator, and to p= —1, in which 
case the surface loads are inversely proportional 
to the distance to the edge of regression, for 
instance in the case of a membrane hyperbolically 
reinforced towards the edge of regression. 

If C and D are zero, in the case p=0, the 
stress flows are zero along the edge of regression, 
and in the case p= —1, the stress flows are finite 
along the edge. The last hypothesis will not be 
examined ; assuming that the surface loads are 
finite in every point, it appears that the stress 
flows along the edge of regression are determined 
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by the factors C and D, i.e.: 








C(A1) 
n,=0; ae ak 2" 
2 1 (12) 
, pC’ (Ai) pC(r1) © D(Ax) 
2=- a 


(A2—A1)? (A2—A)? hemAr 


The generators being tangent to the edge of 
regression, it results that the stress flows nz do 
not appear as forces perpendicular to the edge 
of regression, but as shear flow, the value of 
which is, from Fig. 1, in which the triangle ABC 
is considered as infinitely small: 


sla ye ) 
fijo = —limit 





= —limit 
A2=Al 


A\o—A D(v 
E 2 ~|-- (Aa) (13) 


p p 


If D is different from zero, and C is equal to 
zero, there is a finite shear flow along the edge 
of regression, and more generally, the shear flow 
has the value: 


1 
fiyg = 2—limit |=n:0—a0) | (14) 


A2=A 

2 1 p 
If C is not zero, there is an infinite shear flow 
along the edge of regression, and such a case has 
no practical interest. 


Ill. EDGE CONDITIONS 


From the results established above, the condi- 
tions of statically determined equilibrium of the 
membrane can be deduced. 

Case I. The membrane does not include its 
edge of regression. 
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(IA). The membrane is limited by two 
generators and two lines of curvature (Fig. 2). 

If two conditions are imposed at each genera- 
tor, these two conditions give two relations 
between the unknown functions C(A;) and D(),). 
Thus the arbitrary functions are determined. 
More exactly, if the shear is not given, C(A,) 
appears through a differential equation of first 
degree and it is necessary to give the shear along 
one generator, or a general resultant of shear 
along one of the two limiting curves, thus giving 
the possibility of determining the arbitrary 
constant which appears in C(A,). Moreover, it is 
impossible to impose the value of shear at both 
ends of a generator, because this would give 
two equations for c(A,) and zero equation for 
D(\,). Thus, at the two ends of each generator, 
it is possible to impose one of the following 
groups of conditions: (a) normal flow determined 
at both ends of each generator, and a resultant 


of shear determined along one edge; (b) shear 


and normal flow determined at one end of each 
generator; (c) shear flow at one end of each 
generator and normal flow at the other end. 


Moreover, along a generator which is the edge 
of the membrane, the shear and normal flows 
are determined. Thus a generator which coincides 
with an edge must be reinforced by a member 


accumulating or redistributing the shear flow, 
and also by a beam absorbing the normal flow. 


This beam can be suppressed when there are no 
distributed loads on the membrane, or at least 


in the region of the membrane near the edge. 
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It may also be possible to use a part of the 
membrane along the generator, reinforced if 
needed, to constitute a beam. 

(IB). The membrane is closed, i.e., limited to 
two lines of curvature (Fig. 3). If the shear is 
given at one end with the normal flow given at 
same or opposite and of each generator (Cases 
IAb and IAc), the membrane will be in equi- 
librium. 

If the normal forces are given at both ends 
(Case IAa), it is necessary to verify if the shear 
flow function will be univalued, i.e., if by inte- 
gration of Eqs. (9), m2 will-be the same after a 
complete integration along a closed curve. If 
not, this will mean that the membrane will not 
be in equilibrium under the given loads, and will 
need one or several members along one or several 
generators to absorb the difference of shear flow 
on both sides of the member. If such members 
exist, or if the loads at the end curves are in 
equilibrium with the applied loads, the calcula- 
tion will be made as for the Case IA. 

(IC). The membrane is open or closed, the lines 
of curvature are replaced by two curves which 
do not coincide with the generators. The condi- 
tions in IA and IB are not essentially changed, 
the only modification consists in the fact “that 
the end shear and normal flow not the 
quantities 12 and m2, but combinations of m1, m2, 


are 


and m2. The shear go’ and normal flow mo’ along 








\ 

\ 

\ 
\ 
\ 
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the edge (Fig. 4) will be: 
Ny’ =n, COS* w+ sin? w+2n12 Sin w COS w, 


go’ = (1 —M2) sin w Cos w (15) 


+12(sin? w—cos? w) 


(positive signs as shown on Fig. 4). 

The shear and normal flow can be replaced by 
the flow along the generators and the “general- 
ized”’ shear conjugate to the direction of the edge 
and the generators. Calling these flows n’ and q’ 
we find (Fig. 5): 

No 
n'= =n, COS w Cot w 


sin w ; (16) 
+n sin w+2n12 COS w, 





q’ =qo' +n’ cot w=N12+MN Cot w. 


In this manner it is found that, provided it is 
made use of conjugate flow along the edge, the con- 
jugate shear is not a function of m2; thus it does 
not contain C’ and D, and the results obtained 
under IA and IB remain unchanged, as the 
knowledge of g is equivalent to the knowledge 
of my2 (m, is determined in all cases: n,;= —ZR’). 

On the contrary, if the end conditions are 
given in terms of shear and normal flows q and n, 
a condition involving one of these two quantities 
only does not give the possibility to determine 
C directly. Thus: (a) Conditions introducing 
conjugate flow (n’ and q’) are strictly equivalent 
to conditions along lines of curvature. (b) Con- 
ditions introducing rectangular components (79’ 
and go’) are always equivalent to conditions 
along lines of curvatures in Case IAa (normal 
flow given at both ends). 

C(A,) cannot be calculated directly, and a 
general condition for shear must always be 


a 


Fic. 6. 
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introduced. Moreover, for closed surfaces, it is 
necessary to have a univalued function C(A,). 
It is to be noted that when there are no dis- 
tributed loads, then 


ni=0, n’=n2sinw+2ny.cosw, g’ =m. (17) 
The generalized shear flow is then equal to the 
shear flow along the generators. This greatly 


simplifies the interpretation of the results. It is 
also to be noted that 


Ud 





= MNot+2n12 cot w. (18) 
sin w 

This means that the force m2 to be transmitted 

between two generators is not equal to the applied 

external load corresponding to the distance 

between the same generators. 

(ID). The membrane is limited by a closed 
curve. This curve has two points, at least, where 
it is tangent to two generators at the surface 
(A and B, Fig. 6). In such points, the shear 
may be infinite and it is necessary to check for 
each condition of loading if that is the case. 

Case II. The membrane includes its edge of 
regression (Fig. 7). As it has been seen, C(A,) 
must be equal to zero in order to have a finite 
shear. If along the edge of regression, normal and 
shear flows are zero, and if on the other edges 
(generators or other curves) the shear and 
normal flow are not given, the membrane will be 
in equilibrium and C= D=0 will determine the 
flow in any point. If the shear flow is not deter- 
mined along the edge of regression, D is not 
determined, and C is equal to zero. Thus, in 
order to determine D, it is necessary to have 
conditions of the type: determined normal flow 
along the edge DD’, undetermined shear and 
normal flows along AD and A’D’ (generators). 
Then the membrane is in equilibrium. 
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IV. DISTRIBUTION OF STRESS ALONG GENERA- 
TORS (NO DISTRIBUTED LOAD) 


Equations (10) become in that case : 


C(Ai) 
n,=0, N2= —, 
(A2—A;)? 
’ (19) 
pC’ (Ai) pC(Ai) D(A) 
nN» — + -.-— 


st (A2—A1)? (AeA)? hen 


\, is the are along the edge of regression, while 
Ae—A, is the radius of geodesic curvature r; of 
the curves \, (see Fig. 1). Thus it can be written: 

C(A1) C’(Ax1)  pC(Ax) — D(Aa) 


e+ = — : 
si ri" r,? Yr; 


(20) 


It appears that m2 is given by the relation of 
shear in tapered panels, the apex of the panel 
being at the edge of regression. In order to 
examine m2, consider, instead of the flow, the 
total force acting between two generators sepa- 
rated by an angle da. rida is the distance 
between the two generators at any point. Thus, 
the force is (noting that pda=d),) (see Fig. 8) 


C’(A,) dA, C(A,)dry 
d Fy =nzr,da = D(d,)da+——_—_ + __—_, 











r1 a og 
dC(A1) C(Ai)dAy 
dF, = D(d\)dat+ + —, (21) 
Ti rr 


The factor D(\,) -da corresponds to a direct load 
applied along the membrane between two gener- 
ators. The factor dC(A,):7; is due to the variation 
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of shear from one side of the considered element 
of membrane to the other. The factor C(A,)dA,:r2 
is due to the variation of position of the apex 
along the generators (as compared with a cone, 
in which case the apex has a definite position). 
Assuming that the generators are limited by the 
radii 7;’, 7;’’, with no end loads, then: 
dC(Ai) C(Ai)dA1 
0 = D(A1)da+———-+ es 


Tr) r,’? 


dC(ai) C(Ai) dry 
0= D(A. )da+——-+ — ° 
r\ 








ws 
TY) 2 


Fy, total shear force which is transmitted 
along the generator, is given by 


ery’ 1 1 
Fy2= | nodr, = C(A}) ——-— }: 
VJ r\’ Yr) Yr) 





(22) 
rr" 
C(A1) = Fi — 
7. = F's 
Thus, eliminating D(A,)da and C(A,): 
dC(r;) f 1 1 1 
~-—-— —-—) +Fi.( —+— } =0 
adr, ry’ r;"" ry ry! 
(23) 
dC(A1) rn’ +n 
dd, 7 “=r ‘ 


This relation gives the variation of the distribu- 
tion function C(A;) with the position of the apex. 
The value of D(\,) is 








pF yo . 
D(\1) =— es (24) 
Ke (ry —r1) 
Thus 
dF» Fie ri +r,’ ri’ry’ ; 
nt, 4 1————_+——_|._ (25) 
dr, ry —ry’ r ry? 
The maximum of this quantity is obtained for 
a: 1 
—=—-+ (26) 





Yr) ry’ r;" 


and has the value 


dF; Fr(ri’—rn") Fiuf 1 1 : 
(=) - ‘ ==). (27) 
dy, m Ar,'r"" 4 r,'" : vr,’ 


Calling 








7 
Tr) ° 
——=t=taper ratio; 
ry 
(28) 


L=r,"'—r;' =length of generator, 
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it is found: 





} a i 1 1 (t—1)? 
ry =—; ry =} masts =———.. (29) 
t—1 t-—1 i Hh Lt 
and 
PE we: (30) 
dd, = —dr,|’ =———_.. 
(t—1)? 
Thus 
dt 
(dF), = a Fy.—. (31) 
4 


If the variation in taper is finite, it may be 
found, ¢’ and ¢” being the two tapers, 


VE —J/t 
(AF 2) n= — Fi2— a (32) 
Mtl + Jf/l' 
which, at the limit, gives the relation (31). 
If instead of considering F2, nz is introduced, 
it is found 


noda Fy. Lori’ try! ory’ ry’ 
= |-- + “| (33) 


ry r;* ry" 











dr, ry’ —ry’ 


which is maximum for 


s 3 1 1 1 1 ; 
7 Rage — ahaa ‘ow (34) 
ry lal Yr) i ME so Wied hk 


and has the value 


1 2 1 
("2da) m = -| - (v +—) +5(#+ ) 
Vt 3 t} 


2 1 111 Frye dt 1 
+-(t+--1) |} ————_., 
3 t 9Lt ’ 1 


(12) m= | on (. 4) +2( 04 ) 
912 Jt} 3\ 8 
2(e41) (vet), 
3 t Vi 


n, and F2 are compressions when the shear is 
directed as indicated on Fig. 9. 

When the normal flows m2 at the ends of the 
generator are different from zero [this generally 
happens when the edges do not coincide with 
the lines of curvature, see Eqs. (17), (18) ], it is 
necessary to determine the points of zero normal 
flow (which can be out of the portion of the 
membrane limited by the edges), and apply the 
above formulas between these points. 








(35) 
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e>0 


Fic. 9. 


Thus, being given at both ends of the generator 
(see Fig. 10), 


; n' —2n,2' cos w’ 








te = r ’ 
sin w 
? 7 7 (36) 
n’’ —2n 2’ cos w 
7 
n= ° ” ° 
sin w 


the equations giving C and D become 








dC(A1) C(Ai)dA1 
r'danz’ = Ddat andl ae 
ial at ts 
dc(X1) c(Ay)da 
r'dans"" = Ddat+ + pe 5 
ial re 


Equation (22) remains unchanged while (23) 
becomes 


dC(d;) f 1 1 1 1 
(—-—) +Fa(St—, 
dr, ri ial ri r\ 


(n2'r’ ae no'r'’) 




















= (37) 
p 
and 
ne!'r "2! p 19 pF 2 
Pi) ae 
= ry” 75 it Cas 
Thus 
dF, Fy. ry try ry ry’ 1 
ene yo = sea ais ” ’ 
dy oni’ —nry ry ri” p(ri’’—ri') 


ry’ ne” 
x| nars( 1 ae -) —niiri*(1 --)| (39) 
Y ri 





1 1 \pdF, : 1 1 1 1 
ar “mS ae 
ry Wy dd, ry 1 ri ri 
1 1 1 1 
+nin"(—-— —n2 ry" 5 eat ° (40) 
Lal r) ri ri 
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In order to determine the points of zero value 
of dF,/dd,, the following unknown is taken: 


2/ry—1/r)' —1/r," =x. (41) 
This leads to the equation 
2 woe ot 1 I : 
x* p Fp — 2x(ne YT; Ne; )— pF» —-—;) 
ri ri 
1 1 
+2(nn'+m2"rs")(—-—) =0. (42) 
ri ri 


The values x’ and x” thus found give the 
locations of the points of zero dF.,/d\,. The point 
of maximum dF,/d), is given by x corresponding 
to the average of the two values x’, x’’, and 
corresponds to 














no!" —Neo ri 
xXo= (43) 
om pF i2 
hus 
1 1 nor) —N» ry’ 
aaa 2 i i (44) 
ry OS al pF i2 
For that point, 
dF, 1 1 Fie nol ry’ +n’ ry’ 
sp “Ite ha? 
dr, m ry; r\ 4 2p 
(no"’ry"’ —n2’r)')? 
- (45) 





+ 
| 1 1 

te Fa — ;) 
ry) ri 


So far an isolated generator has been examined. 
The results can be extended to a complete surface 
under determined conditions of loading. Here- 
after, the condition of shear transmission on a 
surface limited by two lines of curvature is 
examined. For that purpose, equations written 
as (19) are used. Repeating the same calculation 
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as before, Eqs. (22), (23) are found and, elimi- 
nating Fy. between those two equations and 
replacing r in function of X, it is found: 


1 dC(aA) ry’ +r,’ ( 1 4 1 ) 
C(A,) dr, 7 rir" “ r)’ ven 
1 1 








= ——_+-___. (46) 
Ao’ Ar Ae” Aa 


By integration, 


" log C(A,) = log (A2’ —A1) (A2”’ — Aj), 
oO 


C(A1) = K (A! — Ai) (A2” — A) = Kr’ 
= Fyory'ry”": (ry —ry), 
C’(A1) + K[(A2’ —A1) + (A2”” —A1) ] =0. 


Then 


(47) 
(48) 


K (2! — 1) (A2”” — 1) 
- (A2—A1)* 
Kry'r," Fyor ry" 
- - , (49) 
ry? r2(ri’—ry') 
p(r2"” +A2’ —2y,)K 
(A2—A1)? 
p(A2’—A1) (Ao —A1) KK D(Ai) 
(N2—Au)? e—At 





ni2>= 





1o=— 








ae [ (Ao! — As) (A2” — Aa) 
(e—da)? 2 2)\A2 2 


oq, D(A) 
= (A2—A1)? |+ eit oe 
Ao—A1 


In order to have n2=0 at the ends, D(\;) = pK. 
Thus: 
pK (ry —ri)(ri’ —r1) 


3 





ni2= 
r\ 


pF y2(r1’ —1r1) (ri — ri) 





(50) 
(ry —ry')ri 
The two formulae (49), (50) thus solve the 
problem of shear transmission. In_ practical 
applications, the problem of shear transmission 
is stated as follows. The normal loads applied 
at the end of the generators being given, deter- 
mine the stress distribution at every point of 
the surface. Let m2/(A1), m2’’(A1) be the values of 
m2 at the end of the generators. Then we may 
write the last equation (19) applied at the two 
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ends of each generator as: 


epC’(A1) = pC(A1) 





(Ax —A 2’ (X = D(r " 
tits M-mv —Ay? se 
a 
C’(A3) C(A1) 
a cc 





Ao’ =A (A2”” —A1)? 
Taking the difference of the equations, 


(A2’ — A1)m2'(A1) —(ro"” —X1)m2""(A1) 


1 1 
= pC’(d)| ——-— — 
Ae — yi Ae —A1. 








1 
+c0)| - :} 
(A2’—A1)?_ (A2”” — A)? 


Writing C(A,) in the form: 
C(A1) = K (Ai) (A2’ —A1) (A2”” — Aa), 


the above equation becomes with the new 
unknown K(A,): 


pd K(X) 


dx 





(A2’’ — Ao") = (Ao! — Ar) 2’ (A1) 
— (o’’ —d1)me""(Aj). (52) 


By integration K(A,) is determined, and _ re- 
placing it in function of C(A,) it is found: 


(A2’ —A1)(A2”” —A1) 
C(A1) =— = re! | fCar-aynzan) 


id ad dd, c 
— (Ao ~Ay)me2’’(A1) F—. = (53) 
p 








Thus 
(A2’ —A1) (A —Aa) 

ee ; 2 Jf Cors—anyms' Ou) 
(A2" — Aa’) (A2—A1)? 


( ” , dX. 4 
— (deo! —A)m2"’(A1) }— (54) 
p 





and by substitution in the expression of m2 and 
determination of D(\;) by one of the equations 
(52), m2 is found. 


(A2”’ —A2)(A2’ — Ax)? 
N»=- No 
(A2’’ — Ao’) (A2—A1)? 
(A2’ — Az) (A2”” — Ax)? 

Ne 
(A2’ —A2””)(A2—A1)? 
(A2’ —A )(A2” —d )p 
— - ~ fCas—aynon) 
(A2—A1)*(A2” — Ao’) 
a7 7 dy 
— (Az —1)M2 (A1) }—. (55) 
p 





“(A1) 





(1) 
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V. PRACTICAL APPLICATIONS TO AIRPLANE 
STRUCTURES 


In actual structures, such as wings or fuselages 
of airplanes, the developable membrane which 
is constituted by a skin-stressed structure is 
divided into segments by redistribution frames. 
Each segment between frames can be analyzed 
by means of the methods developed thus far, 
provided the edge loads are given in such a way 
as to render the segment statically determined. 

From the examination of the edge conditions, 
it appears that it is necessary to determine axial 
load distribution in the direction of the genera- 
tors along both frames limiting a segment. Thus, 
the knowledge of axial distributions at each 
frame renders the structure statically determined. 

Since the distribution of axial loads is not 
known except at the end frames, it is necessary 
to introduce a method of distribution of those 
loads. For structures of elongated shape in the 
general direction of the generators, having 
parallel frames perpendicular to such general 
direction, hereafter called axis of the structure, 
the beam distribution theory is most convenient, 
since it defines axial loads in a section by the 
properties of the structure at that section. 

It is to be borne in mind that such method of 
distribution is arbitrary, at least for two reasons: 
(a) Beam distribution applies only for cylindrical 
beams, and theoretically, for constant shear loads 
only; (b) from the distribution of axial loads 
along a generator which has been established in 
the preceding paragraph, it appears that if a 
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beam distribution exists at two end frames, it 
does not subsist at intermediate sections, which 
shows the arbitrary character of the assumption 
of beam distribution at redistribution frames 
only. The beam distribution at frames is justified 
by its simplicity and by the existence of phe- 
nomena of buckling and plasticity, which prevent 
the determination of effective areas with a degree 
of certainty which would warrant a research of 
actual distributions. 

Omitting the conditions which render a beam 
distribution acceptable in the case of a cylinder, 
it is always possible to check the legitimacy of 
the extension to developable surfaces by calcu- 
lating the axial load distributions at various 
sections between frames. If such distributions 
are sufficiently close to beam distributions, the 
method can be considered as acceptable within 
the degree of accuracy required. If not, the 
method must be the function of 
determining stresses inside a structure 
cially shears 


reduced to 
espe- 
which are compatible with applied 
edge loads. 

The beam distribution method defines compo- 
nents of the loads along the axis of the structure; 
if a is the angle between a generator and the 
axis the quantity thus determined is i=n cos a, 
n being defined on Fig. 5. Substituting in Eq. 
(18), it is found: 


n 


ne= —2n 2 cot w. 


COS @ Sin w 


(56) 


The actual axial flow nz being different from the 
flow 7 defined by the beam theory, there is an 
additional reason of discrepancy between the 
actual distribution and the distribution deter- 
mined by the proposed method. The discrepancy 
can be reduced by multiplying actual areas by 
the factor cos* a sin? w which compensates, in a 
calculation by a method of least work, for the 
ratio between mz and fi and for the actual length 
of generators, but no correction is possible for 
the term containing m2 in Eq. (56). 


VI. GRAPHICAL METHODS OF ANALYSIS 


In the following paragraphs, it is assumed 
that surface loads are negligible. In actual 
problems of stressed-skin structures, the con- 
tinuous distribution of material is, for conveni- 
ence, replaced by concentrated areas; the method 
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space polygon 








force polygon 


Fic. 12. 


is especially accurate when the structure is 
composed of longitudinal stringers (members) 
and a skin, an effective area of which is supposed 
to be concentrated at the members, the skin 
between members transmitting shear only. Then 
the developable panel between two members is 
replaced by a plane panel; thus, the given 
membrane is replaced by a polyhedron composed 
of small tapered panels. The error resulting from 
that substitution consists of the non-coincidence 
between the line of action of shear and the chord 
of the panel; the difference is generally negligible 
but can be taken into consideration. 

The equilibrium of a tapered panel belonging 
to a polyhedron is defined as soon as the shearing 
force along one of its edges is given, when it is 
assumed that no axial force can act along the 
edges. This results from the fact that along the 
edges in the direction of AA’, BB’, there can be 
no load because of the angle between the panels 
(equivalent to m,=0 in the developable mem- 
brane), and from the distribution of axial loads 
which are considered as concentrated loads along 
the lines AA’, BB’ (Fig. 11). Since all loads on 
the panel are pure shearing forces, there are 3 
unknowns (shearing forces on 3 sides), and the 
equilibrium can be obtained by conditions of 
statics. Graphical analysis can be applied to 
that problem in a very simple manner. Replacing 
the panel by a latticed structure by adding one 
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of the two diagonals A B’ or A’B, the equilibrium 
is obtained as follows: 

In A, forces (AB), (AA’), (AB’) must be in 
equilibrium, which defines the values of (AB’) 
and (AB) [see triangle (1) ]. In B’, forces (AB’), 
(A’B’), (BB’) must be in equilibrium, which 
defines the values of (A’B’) and (BB’) [see 
triangle (2) }. For checking purpose, the diagonal 
A’'B can be used, which defines the two other 
triangles of the force polygon. The space polygon 
and the force polygon are similar figures but it 
must be noted that: The side AB of the space poly- 
gon is transmitting a shearing force (AB) which in 
the force polygon corresponds to A'B': To the small 
side corresponds the large force, and reciprocally. 
That property is the result of the well-known 
relation in tapered panels: The shear flow is 
inversely proportional to the square of the length 
of the side. Since the shearing force is equal to 
the product of the shear flow by the length of 
the side, it appears that the force is inversely 
proportional to the length of the side. It is to be 
noted that this is valid for the two parallel sides 
only, and on the non-parallel sides, shearing 
forces are proportional to the sides, while shear 
stresses are constant. Depending on the conven- 
tions of sign for shear and the sense of the end 
forces which produce the shear, the force polygon 
can be placed in two opposite manners, which 
result one from the other by a rotation of 180°. 

When the shearing force transmitted between 
two points along a generator must be determined, 


the graphical construction is the following 
(Fig. 12). 

The equilibrium of the panel being determined, 
the load along A’A”’, for instance, is determined 
as follows: The line A’’B”, parallel to AB, A’B’ 
is drawn and the equilibrium of the panel 
A’B’A"B” is determined in similar manner to 
the equilibrium of the whole panel. The only 
difference is that the starting line is (A’B’). The 
two diagonals A’’B’ and A’B” are drawn and 
the force along AB” is determined (passes 
by a’’b”’). 

If the load between two interior points of a 
generator A’’A’”’, for instance, must be deter- 
mined, the construction is repeated between 
those points, starting from A’’B’’, in the same 
manner as previously made for the pane! 
A'B’A"”B". A final check can be made by 
determining the load on AB, starting from 
A’’B’"”, or reciprocally. Figure 12 is  self- 
explanatory for the sequence of operations made. 

The analysis of a panel is thus defined when 
the shearing force transmitted from an edge 
member is known. That shearing force is the 
sum of the shearing force transmitted by the 
adjacent panel and the difference of the end 
loads on the member between the two panels. 
It is thus possible to proceed panel by panel, 
provided there is an edge panel, i.e., provided 
the membrane is open. If the membrane is 
closed, the same problem of the determination 
of the origin of shear exists as in closed cylinders. 
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The analysis can be made on the membrane 
developed on a plane. However, the use of a 
plane of projection generally simplifies the prob- 
lem and gives more information on the equi- 
librium. Two types of projection can be used. 


- (a) Projection on a Plane Perpendicular to 
the Frames 


The frames appear as parallel lines and the 
axial loads acting on the members are determined 
by a beam distribution if they are not known in 
terms of applied forces on an end frame. It is 
to be noted that if the beam distribution gives 
loads parallel to the axis, the loads thus defined 
can be directly used. The method is the following 
(Fig. 13). The axial components of the loads in 
the members are placed end to end on a line 
parallel to the axis, the forces at one frame being 
placed first, the forces at the other frame being 
placed at the end of the forces at the first frame, 
as shown (A’, B’, C’, D’=forces at frame F’; 
A”, B,C”, D” =forces at frame F’’). 

On panel (1) the acting load is A’+A’’, thus 
the force equilibrium of (1) is drawn between 
the ends of A’ and A”’, as shown. Drawing an 
arbitrary line (a;"a;’) parallel to a’a”, the 
equilibrium of panel (1) is defined as indicated 
previously by drawing through a’, a;”’ parallels 
to the sides and diagonals corresponding to a’’, a’ 
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of panel (1) and closing with a parallel (8;'81’’) 
to the side b’b’’. The inversion of correspondence 
between a’, a,’ and a,’’, a; must be noted. Thus 
the load along b’b’’, in panel (1), is 8,'’81’. Adding 
forces B’ and B”, it appears that the load along 
b’b’’, in panel (2), is 82’’B2’, on the same line as 
B,'’B;’, but determined by the intersection with 
the lines parallel to the frames drawn through 
the ends of B’, B”’. 

From 82’’82', panel (2) can be analyzed in same 
manner as panel (1) and procedure can be con- 
tinued for each panel. 

In Fig. 13, the projections of the loads in 
members and of the shearing forces are, on 
frame F’, for instance, 


member a’a”’ b’b”’ ce’ d'd” 
load axqaey’ BiB.’ y2'¥3' 53’50 
side a’b’ b’c’ cd’ 
shearing force . a8,’ B2'y2' ¥3'5;3' 


Dividing the shearing forces by the lengths of 
the sides, the shear flow along the frame is 
determined. 

This method gives shear flows, but does not 
define the locations of resultants. It has the 
disadvantage of being inaccurate for panels 
nearly perpendicular to the plane of projection. 
Its advantage resides in the fact that it is 
quicker than the following method, and it is 
well adapted for symmetric structures under 
symmetric loads, using the plane of symmetry 
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space polygon 


as plane of projection. It is to be noted that 
when a pure shear equilibrium is examined, the 
force polygon becomes a figure similar to the 
space polygon. 


(b) Projection on a Plane Parallel to the Frames 


In this case components of the end loads 
parallel to the frames are required. A, being the 
projected length of a member n on the frames, 
L the distance between frames, F,’, F,/’ the axial 
loads in the member m at the end frames, the 
projections on the planes of the frames are: 


J. = FA,:L; fa.’ =F,"4,:L. (57) 


Knowing the loads f,’, f,’’, force diagrams can be 
determined by a method similar to that of 
Fig. 13. 

Figure 14 shows the procedure. Loads at the 
ends of a’a’’, called A’, A’, are drawn end to end. 
Their sum is a;’’a;’. The equilibrium of panel (1) 
is determined by drawing through a,’ and aj’ 
parallels to the sides and diagonals corresponding 
to a’ and a’’. The operations from there on are 
similar to those of Fig. 13, but the loads are not 
represented along a base line, and must be 
measured along each member. The directions of 
the loads must be chosen in such a manner that 
in the completed figure the sequence be, for 
instance, D’’, C’’, B’”’, A’, A’, B’, C’, D’, which 
shows that two opposite forces B’, B’’ are meas- 
ured in same direction along their line of support 
8,""3;’. It results that compression forces are 
measured in one direction, while tension forces 
are measured in opposite direction. Thus, in 
Fig. 14, starting from the top, compression forces 
are drawn from F’ toward F”’ and tension forces 
are drawn from F” toward F’. Once the force 
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polygon is completed, the value of the resultant 
of the shears acting along the frames is deter- 
mined. It is the force ado. Its location is deter- 
mined by a classical construction of addition of 
forces, which is represented in Fig. 14, using an 
arbitrary pole R and the forces acting on the 
frame F’. The resultant is shown in position. A 
check can be made by repeating the same 
operation at the other frame. 

By combining various unit solutions, shear 
centers or stress centers can be determined by 
intersection of two such resultants. Then, arms 
can be measured for the calculation of moment 
coefficients for torsion bending and the torsion 
bending constant may be determined. In case 
of closed profile, a cut is always made in order 
to start the calculation, and a unit solution is 
made under pure shear transmission. If in Fig. 14 
the force placed on a’a” is equal to a’a” at a 
certain scale, the space polygon for shear trans- 
mission becomes the force polygon. From that 
force polygon and an arbitrary pole the resulting 
couple can be determined. 


(c) Extension to Non-Parallel Frames 


The equilibrium of a developable membrane 
can be examined between two non-parallel 
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IT-FORCE DIAGRAM UNDER 
SHEAR AT FRAME F , UNIT 














I-SPACE DIAGRAM MOMENT €. SHEAR AT FRAME F' 
FRAME F FRAME F° -.2432 ,7.0270 
2 - : 
{™ ha. | ; Fr | 
ie I ~~ I : | . 




























































































2-3 
04 so sas 
| 
7) -500 3 
| @ SHEAR - 
oor RESULTANT 
' 3-4' 
' 
0405 | 4-5" 
\ = REF. COL (7) 
2230 
TZ -TABLE-SECTION PROPERTIES-LOAD OISTRIBUTION 
0 2) |@)| 4/|6)| © (7) (®) '9) Wo) 
POINT) AREA! z,| Az.) z Az? |Az:'r | Me3 N=-l | TOTAL 
' 4 ©} © |-\.8/1.296/ -.2432/-.7295 | -4000/-1.1295 
2 ' ! 1 |=.8] 064/-.0270/-.0810 | -.1000/- .1810 
3 i 2] .2 2} 004] 0067) .0200/-.1000 |- .0800 
4 A 3/3/42] 144] .0408/ 1215 |-.1000| .o215 
5 3 4/12 | 2.2/1452/) .2230) .6€90/-.3000; .3690 
Oo |1L8/18 2960) 0000 1.0000 
V-FORCE DIAGRAM CORRESPONDING TO 
T-SPACE DIAGRAM SPACE DIAGRAM Iv 
APPLIED LOAD AT POINT P 1.1208 =1910 .-9900 
r + = 
.667 + 335 3690 "Oars 
° T 
I'-2" 
2°*-3' 
3-4" 




















Fic. 17. Equilibrium of a developable membrane in projection on a plane perpendicular to the frames. 
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frames. Special precautions must be taken for 
the distribution of the axial forces along a frame 
which is not parallel to the axis of the surface, 
and the method is simple only in the case when 
such a frame is an end frame, on which applied 
forces are known. 


(1) Projection on a Plane Perpendicular to 
the Frames 


The procedure is essentially the same as 
between two parallel frames, the difference 
between Fig. 15 and Fig. 13 residing in the fact 
that the two load lines giving the values of the 
loads perpendicular to the frames are not parallel. 
Thus, it is necessary to draw two directions of 
sets of lines defining the loads at the two frames, 
respectively parallel to the two frames. The force 
polygon is drawn starting at the intersection of 
the two first lines, defining the origins of loads 
A’ and A”, respectively, and it must end at the 
intersection of the two last lines, which gives a 
control of the equilibrium. The detail of the force 
polygon is the same as in Fig. 13, and forces 
and shears are measured in the same manner. 
It is to be noted that force and space panels are 
not similar figures in this case, although all their 
sides and diagonals are parallel two by two. 


(2) Projection on a Plane Parallel to the 
Intersection of the Two Frames 


The construction Fig. 16 is comparable to 
that of Fig. 14 for parallel frames. It differs by 





the fact that the intersections of each panel with 
both frames, instead of being parallel, intersect 
along the line of intersection of the two frames. 
The construction of the force polygon is identical 
in both cases. As in Fig. 15, space and force 
polygon are not similar. For open structures, 
resultants and points of application can be 
determined in a manner similar to that corre- 
sponding to parallel frames. For closed structures 
an important difference must be emphasized : 
Distribution corresponding to shear without end 
axial loads does not correspond to a torque, but 
to a resultant which is located along the line of 
intersection of the frames. This appears from 
the fact that since end loads are shears, their 
respective resultants must lie in the planes of 
the end frames. As they are in equilibrium, they 
must necessarily be located along the line of 
intersection of their planes. 


(d) Examples of Application 


The basic method represented in Fig. 13 is 
applied in Fig. 17 to a structure composed of 4 
panels and 5 members. The method is worked 
out in two conditions. (1) A shearing force 
applied at an end frame, in equilibrium with the 
opposite force and a couple distributed according 
to the beam theory at the other frame; (2) a 
load applied by means of two members to an 
end of the structure, in equilibrium with reac- 
tions distributed according to the beam theory 
(axial resultant and resulting couple) at the 
other frame. 
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General Solution of Two-Dimensional Problems of Elasticity 


LEON BESKIN 
Consolidated Vultee Aircraft Corporation, Allentown, Pennsylvania 
(Received April 22, 1944) 


HE solution of the general problem of 

elasticity given hereafter applies for a 
two-dimensional simply connected domain with 
an arbitrary boundary and given volume and 
boundary forces. 

Since the general problem for an indefinite 
plane with various nuclei of strain has been 
solved,':? it is possible to find a solution which 
would satisfy the differential equations of equi- 
librium inside the domain without satisfying the 
conditions at the boundary. Then, by super- 
position, it is possible to reduce the problem to 
that of volume forces equal to zero and given 
boundary forces. Let 

X(¢e:, ty, Try) =the stresses in rectangular coor- 
dinates (x,y) which satisfy both volume and 
boundary conditions. 


Zl (ee Cy, r.,) =the stresses satisfying volume 
conditions only. 

X, Y the components of the boundary forces 
per unit length 2(é,, ¢,, 72)) =stresses defined by 
the relation: 


M) 


= T— 2°, (1) 


Those stresses ¥ satisfy the conditions which 
result by taking the difference of the equations 
of equilibrium satisfied by the systems Y and ° 
(see Fig. 1). 

Since = and <° both satisfy the volume 
conditions, stresses > are in equilibrium with 
zero volume forces. 

Boundary forces X, ¥ which are in equilibrium 
with = are given by the relations: 


0 0 
X =X —«a, cos (n, x) — 12, cos (n, y), 


(2) 


= 0 0 
Y = Y—r,, cos (n, x) —o, cos (n, y), 


in which nm denotes the outer normal to the 
boundary. 


1S. Timoshenko, Theory of Elasticity (1934), p. 108. 
2A. E. H. Love, A Treatise on the Mathematical Theory 
of Elasticity (1927), fourth edition, p. 208. 
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Thus, the general solution of two-dimensional 
problems of elasticity is determined by super- 
position of two solutions: One (2°) is that of 
the equilibrium of an indefinite domain with 
given volume forces. The other (=) is that of 
the equilibrium of a bounded domain with given 
boundary forces and zero volume forces. This is 
the problem which is examined hereafter. 

It is known that the stresses » can be expressed 
by means of an Airy stress function ® 34 


ab 


Ox? 


: as (3) 
t-y= - ; 
" axdy 


ab 
¢:;=—; 
oy" 


and that the stress function is a biharmonic 
function: 


VV =0. (4) 


At the boundary, the stress function is com- 
pletely defined by Eqs. (2); calling s the arc 
along the boundary, Eqs. (2) become :4 

















Obdy Phdx _ 
ucllnctpn FY —oa tt 
dy? ds axdy ds . 
. ‘ (S) 
bdx eh dy ne 
dx? ds daxdy ds 7 
n 
7 s 
wa 
Fic. 1. 
3 Reference 1, p. 25. 
4E. Trefftz, ‘“Mathematische Elastizitaetstheorie,”’ in 


Handbuch der Physik Volume V1; Mechanik der Elastischen 
Koerper (1928), p. 109. 
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by integration, introducing the notations P(s), 


Q(s): 
ab d® 
<= f xas=P0) —=— { Yas=-0(): (6) 


Obdy Addx 
- IG i aces ds 
- f (Pdy—Qdx)=f(s). (7) 


The derivative of @ along the outer normal is 
determined by combination: 


db O® dx ——e (2 
dn “a” dy dn 


+P) =g(s). (8) 
=) =€6) 


It can be shown that #(M) at any point M of 
the boundary represents the moment in respect 
to M of all the boundary forces acting between 
a fixed point R and M, while d6(M):dn repre- 
sents the value of the tangential component, at 
M, of the same boundary forces. 

Thus, the problem of two-dimensional elas- 
ticity is reduced to the research of a biharmonic 
function ® satisfying the boundary conditions: 


d®(s) 
?(s)= f(s}; ———=g(s). (9) 
dn 


The function ® can be considered as a potential 
of volume distribution of masses with a harmonic 
continuous density p/2x:> 


Vb =p (10) 
V7u=0. (11) 


Then the general solution for @ can be written 
as the sum of a solution of (10) (u being a general 
harmonic function) and a general harmonic 
function ®,: 


V6, =0. (12) 


The general solution of (12) can be expressed as 
the potential of boundary distribution of masses 
y 8 


b= f y-log r-ds. (13) 
B 


>O. D. Kellogg, Foundations of Potential Theory (1929), 
p. 150. 
5 Reference 5, p. 160. 
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r is the distance between the point where 9, is 
computed and a variable point on the boundary ; 
B indicates an integration along the whole 
boundary. 

At every point, the component F, along a 
direction m of the force deriving from the 
potential ®; is equal to: 


OP, 


F.=—= [ ‘ cos (n, r)ds, (14) 
On Br 


if the point is not on the boundary, and is equal 
to 


v 
F,= | -cos (n, r)ds+xv (15) 
Br 
if the point is on the boundary. 
The general solution of (11) can be expressed 
as the potential of double distribution w :’ 


cos ¢ 
= f=* ds, (16) 
B OF 
in which ¢ is the angle between the normal at a 
point Q of the boundary and the radius r con- 
necting Q to the point K where u is computed 
(Fig. 2). 

Expression (16), valid only if K is not on the 
boundary, becomes 


cos ¢ 
.= [east nw (17) 


if K is on the boundary. 
From relations (10) and (13) the general 


7 Reference 5, p. 166. 
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expression of ® can be written as: 


o= {ff log rudAt flog r-v-ds, (18) 
D Rk 


in which dA is an elemental area of the domain 
inside the boundary and the index D indicates 
the domain of integration. Equation (18) can be 
transformed, using Eq. (17): 


o= ff loc | f * aoe ‘|i 
D B r 


+ [tog r-v-ds. (19) 
B 


It is to be noted that the integration of r-w of 
(17) gives zero, as the boundary is supposed to 
have a volume capacity equal to zero.§ On the 
boundary, the normal component of the gradient 
of is: 


_9® oO ds 
—= ff a" f~ eid 
r 


+f” sn nO 
B 


r 


(20) 


y is the angle between the normal at the point 
where 06:0n is calculated and the radius r. 
Since ® and 0@:0n are known at the boundary 
[Eqs. (9) ], Eqs. (19) and (20) define v and w. 
Equations (19), (20) can be written as follows 
(indicating the variables to be integrated) : 


K)d 
e -ff ” res] f 2° cos 3 (n, QK) 2 hia 
TOK 


+ [v0 log rpedsg, (21) 
B 


OP p 


ff cos (n, PK) 
dnp D pK 
-cos (n, OK )-ds 
x| f Ai “hia 
B TOK 


+f vq" COS (n, PQ) 
B 


TPQ 


dsgtmvp. (22) 


Notation (n, QK) designates the angle between 
the normal at Q and the radius QK. 


® Reference 5, p. 330. 
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In order to determine the unknown functions 
w and » in function of @ and d@:dn which are 
known on the boundary, it is necessary to make 
the interversion of the double and simple inte- 
grals in (21), (22). 

It is necessary to show that such an operation 
is legitimate, i.e., that the singularities of the 
integrands do not make the operation impossible. 

(a) If the point P is out of the domain of 
integration, it is obvious that the interversion is 
possible for domain D’ interior to D, since there 
are no singular points in D’. In the infinitely 
small domain D—D?’ (width h), the integrals 
can be written as follows [omitting log rpx, rpx, 
cos (np, PK) which are all finite ] (see Fig. 3): 


ff aan f 2 cos , OK) 
QK 
~ fas f° cos (n, QI): M2 + no 


chal nnt| fee Oe 


TQH 


In the first transformation, K is replaced by the 
infinitely close point H on the contour, in the 
second the infinitely small distance h and the 
arc of the boundary S are taken out of the 
integral and some of the terms are replaced by 
their maxima. It is known that 


p= cos (n, , QH)dsq 


-dsq 


(23) 


(24) 
vee 


is bounded if we is bounded. Thus the factor 
multiplying / is bounded and the product tends 
toward zero. 





Fic. 3 
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CS 





If the integral is interversed, the point Q is a 
fixed point in the first integration, which is made 
in respect to K. The integration outside a semi- 
circular infinitely small domain surrounding Q 
has no singular points, and, inside that element 
the first integration is equivalent to: 


sis dg-r-dr 
we — —-—-—— 
ae g:dg-dr=wa:2r (25) 


and tends toward zero with 7; thus there is no 
singularity for the integration in respect to Q. 
(b) If P is inside the domain of integration, 
the results under (a) being established, in addi- 
tion it is necessary to verify that the integrals: 


cos (n, PK)- ‘dAx 
ff tow renda x: [f- ie (26) 
PK 


(with A interior to the domain), have a signifi- 
cance. This result is a consequence of the fact 
that in the infinitely small domain surrounding 


ff log rdg-r-dr=2xf log r-r-dr, (27) 
0 
dg-cos g:r-dr 
(p=, - 


both quantities being infinitely small with r. 

(c) If Q is on the boundary, the first method 
of integration has a significance, from results (b), 
since the simple integral has a significance. 

The second method of integration must be 
examined. This reduces to the study of the 
integrals : 


cos (n, PK) 
h= f oaisf f log rex—— (29) 
PK 


I= foots f f wh hos sinks 
u 


’PK*'TQK (30) 








when Q is close to P (the simple integral being 
effected along a small arc y), since, when Q is 
distant from P, there is no difficulty, as integra- 
tions around Q and P, respectively, can be 
effected. 

When Q tends toward P, it is sufficient to 
consider a domain limited to a circle of radius R 
around P (see Fig. 4). 
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Let rex=p, it appears that integral J, (29) is 
finite when Q comes in P since: 


l 
Ii<wm fas ff ~Fodvde 
= wm f ds-x- J log pd p, 


(31) 
I, <wm-sm4(R log R—R). 


It is possible to determine an upper limit for 
integral J, (30) when Q tends toward P. 

The domain of integration (semicircle having 
P for center and MM’'=2R for diameter) 
divided into several domains having common 
parts and the integration is made separately in 
each domain. 

It is obvious that 


(MTM' al 





<(1)+(2) (32) 


in which 
(1)\=(MEJIUNM); (2)=(ETIJE). (32a) 


In domain (1), p=rpx<rex and the double 
integral is smaller than: 


Sf Sf os 
a) & 


Replacing domain (1) by (MEJIN’UNM) of 
greater area and integrating: 


dAx 2 
ff — <a R—1)—y(log y—1) ]. (34) 
(1) < 


'QK >-, (35) 
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thus 


dAx 2dAx 
ftp fe 0 
(2) (2) ¥ 
Replacing the domain (2) by the semicircle 
EQIJE: 
dAx 2 ff” 
ff —— <— dp=2r. (37) 
() Vv v9 


Thus, the total integral in domains (1) and (2) 
is smaller than: 
2x 
2x+ : [Rilog R-—1)—y(log y—1) ]. (38) 
That expression is finite 
respect to y. 


when integrated in 
Thus, the second method of inte- 
gration has a significance, and the interversion 
of the integrations is possible. Then, the equa- 
tions at the boundary can be written as follows: 


br= f wo: dsg 
{ff cos (n, QK)-log rpx: ‘dA ‘| 
QK 


ia rpg@dsg, (39) 
B 








OPp 


aus = f w0-d50 
onp B 


Sf cos (n, QK)-cos (n, PK)dAx “| 
TQoK*TPQ 


+f: vg: cos (n, PQ)dsq 


pe 





+rvpp. (40) 
Double integrals in (39), (40) are functions of 
the domain of integration and of the points 
(P, Q) and can be determined for any couple of 
points (P,Q). They will be designated by 
Q(P, Q) and ¥(P, Q); it is to be noted that the 
kernel 2 is symmetrical in (P,Q) while the 
kernal V is not symmetrical. Equations (39), (40) 
can be written, replacing P by V in the second 
equation : 


' tr= f vlog req‘ dsq 
B 
+ [wo 0P, Q)dsg, (41a) 
B 
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O® cos (n, VQ) 
( ~) = | v0: — Pisin 
Ons y B rve 


+ f oo¥( V,Q)dsq. (41b) 


The problem is thus reduced to a system of two 
integral equations. The solution is obtained 
multiplying (41a) by an arbitrary function 
M(P, V), integrating in respect to P, and adding 
the result to (41b). 

Then, provided that 


fel wee, V) f oon, Q)atse| 


+ f wo¥(V, Q)dse=0 (42) 


the system (41) is transformed into 


Ob 
( ) + [emer V)dsp 
Ons y B 
= fase] er, V): f re log rvats| 
B B 


VQ) 
+ fer aset en. 
B Tv@Q 


(43) 


Equation (42) can be written 


feeaisd wv, 
E 


Since it is an identity in V, Q for any function w, 


M(P, V)Q(P, 1 =(). 
+f re O)ds»| 0. 


it is equivalent to 


(V, o+ fur, V)Q(P, Q)dsp=0. (45) 
B 


Integral equation (45) defines a function M(P, V) 
in a given domain D. The function M(P, V) 
being determined, Eq. (43) can be written as 


O® 
i) &p-M(P, V)-dsp+ (=) 
B on V 
= frase f log reg: M(P, V)dsp 
B B 


Phos (n, VQ) 


rve@ 


= |. (46) 
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Let 


f log reg: M(P, V)dsp 
B 
cos (ny, VQ) 
4+ V 


Tve@ 





=A(V,Q). (47) 


Equation (46) becomes 


ob 
fre -M(P, V)-dsp+ ( ~) 
Ons y 
= franc, Q)-dsgtmvg. (48) 


Equation (48) is an integral equation which 
determines v. Any of the Eqs. (41a), (41b) 
determines w when v is known. Thus the solution 
of the general problem of elasticity in two 
dimensions is reduced to the solution of several 
integral equations. The method can be sum- 
marized as follows. 


PROCEDURE FOR THE SOLUTION OF TWO DI- 
MENSIONAL PROBLEMS OF ELASTICITY IN A 
GIVEN DOMAIN (NO VOLUME FORCES) 


(a) Calculations to be made in the domain 
independently of the values of the applied 
boundary forces. 

1. Calculate the integrals: 


cos (n, OK )-log rpx-dAx 
Q(P, o-ff — e eee - ni (49) 
D 


v(P, o-ff cos (n, PK) -cos (n-QK) -dAx 
D 


’PK’TQK (50) 





2. Solve the integral equation for the function 
M on the boundary: 


v(V, a+ f M(P, V)Q(P, O)dsp=0. (51) 
B 
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3. Calculate the integral: 


Vi 
A( V, Q) = Ler onc 


log rea 
TVQ . 
-M(P, V)dsp. (52) 


(b) Calculations to be made along the bound- 
ary for each loading condition, i.e., for each 
group of values of ® and 0@:dn along the 
boundary. 

1. Solve the integral equation in v: 


OP 
[ ee ue, Q)-dsp+ (=) 
B On] vy 


- f vq: A(V, Q)dsqo+nvg. (53) 
B 


2. Solve one of the two following integral 
equations in w: 


br f vq: log rpq'dsq 
B 


4 f wo 0P, Q)-dsg, (54a) 
B 


(=) f vg:cos (n, PQ)-dse 
On P B TPQ 


— r= [ wo-W(P, Q)dsq. (54b) 
B 





3. Calculate the potential » inside the domain 


D: 


we Bcd. Ataf (55 


TQK 


f wae'cos (n, QK)-dse 
KK>= 
B 


4. The stress function ® is defined by the 
relation 


off log rik: ux:dAx 
D 


+f VQ log rrq‘dsq. (56) 
B 








On Van der Pol’s and Related Non-Linear Differential Equations 


J. SHOHAT 
University of Pennsylvania, Philadelphia, Pennsylvania 
(Received April 15, 1944) 


HE present paper is closely related to a 

previous one by the author.' First, we 
correct a consequential error in the said paper. 
Secondly, by a simple application of Parseval’s 
formula we are able to establish directly the 
following fundamental result. If the equation 
(d°u/dt?) —eF(u)(du/dt)+u=0—generalization 
of Van der Pol’s equation—has periodic solu- 
tions, of period close to 2x, then, under certain 
general conditions concerning F(u), such a solu- 
tion u(t) is strongly dominated by the first term 
of its Fourier expansion, that is, not only u(t), 
but also u’(t), u’’(t), u’’’(t) can be approximated 
closely, for all ¢, by the said first term: 
harmonic oscillation- 


simple 
and its corresponding de- 
rivatives. An algebraic equation is set up for the 
determination of the amplitude A, of this 
harmonic oscillation, and Van der Pol’s result 
(A,~2) is confirmed. Thirdly, for Van der Pol’s 
equation we set up expansions of u(t) and its 
frequency v in power series, in terms of ¢,0<e<1. 
Finally, for the same u(t) and v we obtain, for 
the first time, we believe, other series expansions 
in terms of e¢, for all values of ¢, large and small. 
The values of v computed from the first few 
terms of the new series agree quite well with 
experimental results. The method employed is 
applicable to more general equations. 

Corrigenda. |n my paper referred to above the 
statement on page 43: “--- eP(t)|, --- in (31) 
and (31.1) are very small for all ‘=0”’ is incorrect. 
In fact, even a small change in the frequency 
of a harmonic oscillation may affect the oscil- 
lation considerably as time goes by. This simple 
fact had been overlooked by the author. Thus, 
the conclusions drawn on the basis of this state- 
ment, and similar ones applied to the non- 
homogeneous equation, do not hold. But what 
about the relation u?(t)+u(t)~4 [Eq. (36) ], 
which holds according to Van der Pol and which 
was confirmed on the differential analyzer? We 
show their validity below by a purely analytical 
treatment. 


1 J. Shohat, J. App. Phys. 14, 40 (1943). 
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A. PARSEVAL FORMULA 


1. Let f(f) be periodic, of period 7, admitting 
a Fourier expansion 


; ais . 2n 
f()~—+E a, cos nvt+6, sin nvt, v= ~ 
n=l 
(1) 
a, 2 


al cos ‘ 
| f(x) . nvxdx; n=0,1,2,°->. 
0 


a F sin 
Assuming the existence of /o7f?(x)dx, we have, 
whether (1) converges or diverges, the classical 
Parseval’s formula 


2 pT ag * 
=| f?(x)\dx=—+}D a,?+6,’. (2) 
r: 0 2 n=l 


This formula—a very powerful tool in the theory 
and applications of Fourier series—ts the basis 


for the subsequent discussion. Consider the differ- 


ential equation 


dues du 
ru Fees, e>0. (3) 
at 


We make here the following assumptions: (7) 
F(u)- is continuous for all finite «; (i) F’(u) 
exists; (iii) F(u), F’(u), hence, also fo" F(u)du, 
So'uF(u)du, are bounded if u remains in an 
arbitrarily fixed finite interval [all these condi- 
tions are fulfilled if F(u) is a polonomial ]; (7V) 
under certain initial conditions the differential 
equation (3) has a continuous periodic solution, 
of period 7, with frequency 

2r 


y=—.- (4) 
T 


We then have the Fourier expansion 


ao = i - 
u(t)=—+> a, cos nvt+b, sin nvt, (5) 
n=1 


2 Sufficiently broad conditions for the existence (and 
uniqueness) of periodic solutions for differential equations 
of the type (d’u/dt)+f(u, du/dt)(du/dt)+g(u)=0 have 
been given in a recent paper by N. Levinson and C. K. 
Smith, Duke Math. J. 9, 382-403 (1943). 
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2 rT cos 
-—{ u(x) . nvxdx; n=0,1,2,---. (6) 
T Jo sin 


Qn 
bn 


The expansion (5) converges for all ¢, since u(t) 
and u’(t) (see my paper referred to in reference 1), 
hence, also u’’(t), u’’(t) [by virtue of (3), where 
F(u) satisfies the conditions stated], are con- 
tinuous and bounded for all ¢. 

Integrating (3) from 0 to T gives: 


|T 


T oT 
u'(t)| -ef Fudu+ | udt=0, 
10 0 0 


‘ T u 
(w—eg(u) e+ | udt=0, (ows f F(udu), 
0 0 


whence, since u is of period 7, 


T 
{ udt = (0).* 
eo 


Thus, in the expansion (5) a9=0, so that (5) 
becomes 


u(t)=>- A, sin (nvi+a,), An=(a,.2+5,2)'. (7) 
1 


r= 


Furthermore, 


u'(t)=>° A,nv cos (nvt+a,), (8) 
n=l 
u(t)=—> A,n*v* sin (nvt+a,). (9) 
n=l 


(These expansions also converge for all t, by the 
above reasoning.) Multiply (3) by u, also by 2’, 
and integrate from 0 to 7. We get: 


” T T 
[unu’—ep(u) |o — [ watt [i udt 
Jo 


=0(4u)= | uF(u)du), 
0 
T T T 
[ 2u'wdt—26 f F(u)u"dt f 2uw'ar=o, 
v0 0 0 
aT oT 


| wrdt= | urdt, (10) 
e/ 0 


7 
f F(u)u"dt=0 (€¥0). (11) 





* For the more general equation (d?u/dt*) — «F(u)(du/dt) 
+g(u)=0 we get similarly /o7g(u)du=0, as a necessary 
condition for the existence of periodic solutions, so that 
g(u) must necessarily change sign within a period. 
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The relation (10) is fundamental. 


Thus, F(u) necessarily changes sign within a 
period ; hence, (3) cannot have periodic solutions 
if F(u) keeps a constant sign for all «; moreover, 
if uw; is the smallest numerical value of u, where 
F(u) changes sign, then necessarily 
“maximal amplitude” of w=max.|u(t)| > | a1}. 

O<t<T 
We now turn to (10). Apply Parseval’s formula 
(2) to either side, making use of (7) and (8): 


(12) 


> A,2(1—nv) =0. 

n=1 
If A,*(1—v’) <0, then A,?(1—n?y*) <0, for n=2, 
and (12) is impossible. If A?(1—v?)=0, then 
A,2(1—n’v)=0, n2=2, and (12) requires 
A,2(1—n’?v*?)=0, for n=2, which, barring the 
trivial solution u(t) =0, leads to the uninteresting 
case A,=A;3=:--=0, u(t)=Aj, sin (vt+a), 
which correspond to e=0 in (3). It follows that 
if e¢~0 in (3), the first term in (12) is necessarily 
positive, so that 


v<i, T>2r. (13) 


2. O<e<I/. For e=0, all solutions of (3) have 
period 27, i.e., v=1. Assume that for 0<e«<1, 
the period is close to 27, that is, y~1 (this holds 
under certain general conditions). How does the 
periodic solution behave? This question is of 
particular significance, since it was shown (Van 
der Pol, Liénard, Levinson and Smith) that, 
under certain general conditions, the periodic 
solution is unique (disregarding a time-shift) 
and that all other solutions approach it asymp- 
totically, as t+ «. Take e so small that 


v2 
v=v(e) >—. (14)* 
2 
Then 
n’y?—1>0, for n=2, (15) 


and (12) yields 


A?(1—»*)=A27(4v—1)+A3*(9r?—1)+---, (16) 

*See below, Part B,. where we treat Van der Pol’s 
equation by means of power series in ¢ and derive an 
estimate of v(e). 
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where all terms are positive, whence 


Ai(1 io — . 
n< ez 2, (17) 
(n*y? — 1)" 
which shows the magnitude of the coefficients 
A», As, «+: in the Fourier expansion (7), com- 
pared to A,. More precisely, write 


u(t)=A, sin (nvt+a,)+OQ(1, €), 


D2 


Ot, eo) => A, sin (nvt+a,) 
siti sin (mvt+a,) 
=> A, (n*v? — 1). 
n= (n%v? — 1)3 
Apply Cauchy inequality ((2a,b,)?=2a,72b,") 
to the right-most term above. This gives: 


z£ x 1 } 
| Q(t, €) | 3(x A,2(n*v—1)- —— ) 
n=2 nae ny? — 1 


for all ¢, and since, by (14), 
1 
Pm— > —2>(s— i}, ea 2, 
r 


we get finally, 
that 


making use of (16) and recalling 


 - 


— nai 


awl 9° 6 


—— (vt+ai)+Q(t,€), (18) 


| Q(t, €) | <a) 2ay (1—,*)} 
“ =0((1—»*)!).* 


This holds for all ¢, under assumption (14). 
Rewrite (18) as 


u(t)=A sin y+2A,(1—vr*)!-0(t, €), (y=#+ay)) 
| A(t, €)| <1 for all ¢. ) 
u(t) thus differs but little from the simple 


harmonic oscillation A, sin (vt+a;) for all ¢. 
What about its derivatives? 
3. By means of (19) and (2) we get readily: 


She -f wdt= fw 


TA 
=—(tm), 


(20) 


m=O((1—»*)!), 


° * The now customary notation x =O(a@) means: x is of 
the same order of magnitude as a; for example, (2n*—3)-} 
=O(n-'), as n> ~, 
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> A,2n*v?=O0((1—v?)!). 


n=2 


[It will be shown below, 
that A i™ 2. ] 


(21) 


on the basis of (18), 


> A,Z=Arm= 


n=2 


O((1—v*)#). (22) 


Again from (3) 


r T oT 
f w'rdt—ef F( u)u'u'at+ | uu"dt=0. (23) 
de 0 0 


By (2), (9), and (19), 


[= f wdt= 5 A,?n'y'4, (24) 
n=1 
u(t)| SA,[1+2(1—v*)!]<3A,, for all ¢. (25) 
Now, by assumption, (25) implies | F(u)|< a 
certain constant C, for all ¢. Apply Schwarz 
inequality 


b 2 eb b 
(( [[ sorpterae) =f fie(x)dx- [ f2(s)dx) 


to the second and third integral in (23). This 
gives, by (22): 


if “Fawu'at| =c( fr wat wd it) 


<C[A?-T/2-(1+m)]}}-T, 


| pT l T T j 
f ww"at =( f u*dt - f w'"dt) 
0 eo 9 


S[(A(T/2)(1+m)]!- 7. 
Equation (23) now becomes 
I—[A;?-T/2-(1 +1) }}(—€6,C+ 62) -['=0, 


6;,2=const., | 6:2) =1, whence 


¥ A2n'yt=A 2(1—v'+0(8))=O0(8), (26) 


5=e or (1—»*)!, whichever is of smaller order of 
magnitude. Thus, As, A3, ---, are even smaller 
than given by (21). By virtue of (26), we now 
get: 

u'(t) = 


Aw cos (vt-+a;)+Q'(t, €) (27) 


Q'(t, €) at A,nv cos (nvt+a,) = O(6') for all ¢. 


n=2 
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The estimate of Q’(t, €) is derived by a reasoning 
quite similar to that employed above in esti- 
mating Q(t, e). Combining (18) and (27), we get: 


u>(t)+u’*(t)~A,’, for all t, error =O(6'). (28) 


The approximations to w’’(t), u’’’(t) are now 
furnished by the differential equation (3) itself, 
namely : 


u’’(t)= —Ayw* sin (vt+ai)+Q''(t, 6), 
| "(t, | =O(6 
QO''(t, €) ) (29) 
u’’’(t)= —A* cos (vi+a,)+Q""'(t, €), for all ¢. 
1Q’'"(t, €)| K1. 


Thus, if «1, the periodic solution u(t) of the 
differential equation (3) and its first three deriva- 
tives can be approximated, for all t, by the simple 
harmonic oscillation A, sin (vt+a,), which repre- 
sents the first term in the Fourier expansion of 
u(t), and its corresponding derivatives. 

4. We proceed to evaluate A,. Here we make 
use of (11). Assume F(z) is a polynomial: 





















F(u) =aptawu+-:+-+a,u”. 





(30) 





By (18) and (27), for any fixed positive integer 
s and for all ¢: 





us=A; sin’ y+ns, = n.(t, €-)=O((1—»*))), 


u”=A ; cos? y+0(5}). 






Substitute this into (11): 





T 
a,A ;**? sin* y cos? ydt+0(6}). 


0 


0=> (31) 


But 





eT 
| sin?'+! y cos? ydt 
0 





T 
=— { (1—cos? y)' cos? yd(cos y) =0. 
o7 0 





We conclude that odd powers of u in F(u) have 
no influence on Ay, within our approximation, 
so that we may keep in F(u) even powers only. 
Returning to (31), where we may set aj=az3 
=---+=0, and using the simple formula 



































1-3---(3~1) 


aT _ 
| sin*' y cos? ydt = ———-- ‘=i, 


0 2i+1 2-4-6---27 





VOLUME 15, JULY, 1944 








we get the desired equation for evaluating A,: 


‘Ss 1-3---(2i-1) zi 
ao 2; i+1 





=0, s=A r. (32) 


imt 2:4-6---2% 
Note that in (32) the coefficients are independent 
of « and 7. The same equation was obtained in 
an entirely different way by Liénard.? We 
conclude that A,;—the square of the amplitude of 
the fundamental harmonic of u(t)—ts close to a 
root of Eq. (32). Hence, a necessary condition for 
the existence of periodic solutions for Eq. (3), [F(u) 
having the above properties | is: (32) possesses at 
least one non-negative root. 

As to the phase a; in (18), it is of no special 
significance, for we can replace ¢ in u(t) by 
t+arbitrary constant, since time does not enter 
explicitly in (3). We may set, if we wish, a,=0. 

4. Illustration. F(u)=1+au—u". Here (3) 
possesses a unique periodic solution (cf. the paper 
quoted in reference 2). Equation (32) becomes 


1-3---(2m—1) 2 


2-4-+-2m m+1_ 





whence 
2-4-6---2m 


1/2m 
Ay~| ——_——————_-- (m +1 
, Feserey 89) 


The cases m=1, 2 yield, respectively, A1~2, 
A,~8'. Both results agree with those previously 
derived by the author* (a=0, m=1 corresponds 
to Van der Pol’s case). 

Remark. From the classical results of Bendix- 
son‘ it follows that, for «>0, any solution of 
Van der Pol’s equation approaches indefinitely, 
for i, a stationary state represented by its 
unique periodic solution, which, in turn, if 
0<e<1, can be closely approximated, for all ¢, 
by the simple harmonic oscillation 2 sin (vt+ a), 
as shown above. This is the fundamental result 
of Van der Pol. Thus, our periodic solution is 
stable. Its behavior and analytic representation 
are of fundamental! importance for all solutions. 


3A. Liénard, Rev. Gen. de l’Elec. 23, 901-946 (1928); 
Cf. also: N. Kryloff and N. Bogoliouboff, Introduction to 
Non-Linear Mechanics, free translation by S. Lefschetz 
(Princeton University Press, 1943). 

* The estimates (18) and (27) which hold for all ¢ 
justify the method employed by the author to evaluate 
A, in his previous paper. 

4L. Bendixson, Acta Math. 24, 4-88 (1901). Cf. also 
Levinson and Smith, and Liénard, references 2 and 3. 
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B. 0<e<1. APPLICATION OF POWER SERIES 


We exhibit the method, essentially due to 
Lindstedt,’ on Van der Pol’s equation. 


d*u du 
—e(1—u*?)—+u=0. 


ra (34) 
df* dt 


It applies to the more general equation (3), 
F(u) being a polynomial. 

5. We start with the expansion of the periodic 
solution u(t) and of the frequency »=2x/T in 
power series of the parameter ¢. However, a 
straightforward substitution into (34) of a power 
series in € for u(t) is of no avail for it brings out 


nt. We 


cos 


so called ‘“‘secular terms,” like ¢t-. 
sin 


proceed as follows. Introduce a new independent 
variable 


x=2nt/T, (35) 


where 7 is the original period of u(t), so that 
becomes a function of x, of period 2x. Equation 


(35) transforms (34) into 
ru’ —ev(1—u?)u’+u=0, (36) 


where primes denote differentiation with respect 
to x. Here we substitute 


v=yotrvetry2e+:-:, 
u=Po(x) +Wil(xdet+yo(x)e+ et 


(37) 
(38) 


Yi=const., 


a eu = Welle) tdi (et del (xe+- oe 

dx 
where the y’s are periodic functions, of period 27, 
which we take as represented by Fourier series. 
Since ¢ does not enter explicitly into (34), we 
may fix the ¢ origin, ‘=0, so that 


du(t) 


=() for +=0, 
dt 


u(O)=max u(t), 
OStsz. (99) 


Then in (36) u’(x)=0 and u(x) is maximum for 
x=. Letting 


u(O)=aptaetarce+:--, (40) 


it results: 


Vi(O)=a;, W/(0)=0; i=0,1,2,---. (41) 


5 Lindstedt, Mém. Acad. Sci. St.-Pétersbourg 31 (1883) ; 
cf. also Poincaré, Les méthodes nouvelles de la Mécanique 
Céleste (Paris, 1893)? Vol. 11, where Linstedt’s series, 
“interesting and useful,”’ are elaborated in great detail. 
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Since (34) does not change if we replace t by —1 
and « by —e, hence, x by —x and ¢ by —e in 
(36). This suggests the determination of the y—s 
so as to have, for 1=0, 


¥2i(—x) =P2(x), Prigi(—x) = — Pau 
eoi(x) = >>, a; cos lx, 
Poi4i(x) = >0, by, sin lx, 
> lb =0, 


provided this leads to v(—e)=v(e) [see formula 
(55.1) below ]. Then a) coe EY = 73 °°° 
=0, 


i(X), 2) 


(43) 


u(0)=ap+tarxe+aye'+---, 
v=yotyet+yse't+:::. 


[Part C obviates the necessity of showing that 
all Y—s can actually be determined in this way. | 
This greatly facilitates our computations. 

We now substitute (38) and (44) into (36). 
Annulling various powers of ¢€ yields a sequence 
of differential equations for the successive de- 
termination of Yo, ¥1, v2, ---, all of period 27, 
subject to the initial conditions 


W2:(0) =0, 


(44) 


Wo;.1(0)=0; i=1,2,3,--- (45) 


[ the remaining initial conditions (41) are satisfied 
automatically by (44) ]. Here are the differential 
equations for Wo, ¥1, v2, v3: 
yoo +Yo=0, 

Vi + =(1 — Wo" )Wo', 
Yo’ +ye2 — yao’ + ( 1 — Per)’ Pit 2Woo'V1, 
v3 +3 == Qy oi’ + yo 1 — Po? )Wo’ = QWob wi’ 

+ (1 — Yo?) Po’ — (W2+ 2PoWo) Wn’. 


(46) 
(47) 
(48) 


(49) 


The left-hand member in each of these equations 
is of the form y’’+y=0, whose general solution 
is A cosx+B sin x. Thus, the constants a2;, Y2i 
will be determined by the disappearance from the 
right-hand members in Eqs. (46)—(49) of terms of 


cos 
the type const. x *? for the presence of such 
sil 


terms brings out secular terms. This is the main 
feature of Lindstedt’s method. We get first 


Vo= 1, Yo= ao COS X, 


(50) 


y+ vee aoe = : ao* 3x 
1 tyhi= ; ao sin +— sin 3x. 
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Hence, 

ay* , 

——ay=0, ap=2 (for aa~tmax), 

Wo(x) =2 cos x. (51) 
Equation (50) now becomes 
Wi’ (x) +(x) =2 sin 3x, 

and vields, by virtue of (44), 

¥i(x) = 3, sinx—!, sin 3x=sin?x. (52) 
A somewhat simpler method for determining yy; 
and other y's is as follows: 


Integrate (47) from 0 to x, writing yi(x) 
=, b, sin lx, with , /b,=0. This gives: 


1 b, 
pw b(1- +) cos lx+>)> — 
; l il 


8 cos* x\ |" 
={ 2 cos x— “— = 24— 24 cos 3x, 
‘ 10 


1=1,2,3; bo=0, b=—4, 0:=—3b3= 34. 
The same method, applied to (48), yields for v2 
and y2(x) =, a, cos Ix: 
v2= — Meo (53) 
Yo(x) =A cos x+ 3;6 cos 3xn— 5956 cos 5x. (54) 
In order to find A and ¥3(x), apply the same 
method to (49), rewritten as 
v3 +3 = — Zvi’ + 72(Wo— Wo?/ 3)’ 
+L(1—Yo?)¥2—pow? )’. 
We find A = —!g. Thus, 
u(x) =2 cos x+esin® x 
+e (— 1¢ cos x+ 3{¢ cos 3xn— 596 cos 5x) +---, 
2+ Vle+ +++ SUmax=U(0) =2+ o6@+---, (55) 
2 


T=—=2r(1+ce+---). 


Vv 


pu l—Kieet+---, 


[The left-hand inequality for tax follows 
from (25). ] We may check our computation by 
the following relation derived from (36), (51), 
and (54): 


u(0) 2+(A+!366)P+--- 


u'"(0) 24+(A+3%e)e+:-* 


We conclude that 


(1—v*)*=O(e), (56) 
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and this enables us to make more precise the 
estimates for u(t), u’(t), ---, arrived at in Part A. 
Namely, we can replace in formulae (20), (21), 
etc., O((1—v)*), O(6), O(64) by O(e), Oe), O(e*), 
respectively. 

6. The case €>1. Many of the foregoing results 
hold, like formulae (7)-(11). Thus, we again 
have (12), with the resulting inequality 


v<i, T>2z. 


However, here v becomes very small, as e€ 
becomes very large. In fact, Van der Pol and 
Liénard have shown that 


T~1.6:--€, ve~3.89--- (e—-2). (57) 


It follows that in (12) several terms at the 
beginning are positive. Assume, for example, 
1 1 
———_ <FS—. 
(k+1)? k? 


Then we write, in place of (16), 


A(1—v*?)+A2(1—4r*)+---+A2(1—k*v?) 





oe 
= > A,?(n’*v—1), 
n=k+1 
where all terms are non-negative. This suggests, 
as was stated by Van der Pol, that in the 
Fourier expansion (7) for u(t) several terms from 
the beginning dominate (which may account for 
the complicated graphical picture, as given by 
Van der Pol), and the analysis of Part A becomes 
greatly complicated. The power series expan- 
sions, as developed above, obviously cannot be 
applied here. It would be natural to attempt 
expansions in negative powers of e, but the very 
first equation for ~o(x) [see (46)—(49) ] yields a 
non-periodic complicated function, and this bars 
the determination of yj, we, ‘, as periodic 
functions. This difficulty has been overcome in 
the last part, C, by means of a new type of 
series for u(x) and ». 


C. SERIES EXPANSIONS FOR u AND yv, 
FOR ALL VALUES OF e 


7. The method which follows is applicable to 
any equation of the type (3), F(u) being a 
polynomial, where (Liénard) v~1/e for large e, 
[see (57) ]. We apply it in detail to Van der Pol’s 
equation (34), it being the most interesting 
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equation of this type, which played an important 
part in applied sciences, in developing the theory 
of relaxation oscillations, and indeed, the theory 
of non-linear differential equations of physics 
and engineering. 

Again we transform (34) into (36) by means 
of the substitution (35). The main feature of 


our method is to consider the equation 
vu’’ —vep(1—x*)u’+u=0, (58) 


where p is a new parameter in terms of which we 
expand u(x) and v; ultimately p is replaced by 
unity. As before, we set 


u(x)=>> vi(x)p', (59) 
y; periodic, of period 27, 
u'(0)=0, u(O)=Unax; 
y./(0)=0, 71=0,1,2,--: 
Guided by (57), we now set 
1 
eng 
v=(1+h—p)(yotvipt+y20"°+---), (60) 


where h and y; are constants. We note that p=1 
implies e=1/h, «1 for h very large, €>1 for h 
very small, and we may anticipate that our 
series will serve for all values of ¢€, large and small. 

Now substitute (59) and (60) into (58), annul 
various powers of p, and obtain a system of 
differential equations for Wo, wi, --- [see (46) 
(49) |, the first of which is 


Rey eho'’ +o == (), k = 1 +h. (61 ) 
We get, since Wo(x) is of period 27 and yo’ (0) =0, 
1 1 


Jeo ~= 


—-—, (62) 
k 1i+h 


Wo = ay COS X. 


The equation for ¥; now becomes: 
Wi +h s— 2(ky1 —vo0)Wo’ + rol 1 — Po") Wo’. 


As in Part B, we seek to avoid secular terms, and 
the calculations are quite similar. Setting p=1 


yields: 
~ 20 
G3) 
€ 


(x) =— 
v1 


+] € 
SPR ea And you"; 


€ 1+e 


sin® x. 
€ 


Yo(x) =2 cos x, 





The equation for y2 now becomes: 
po” +y2 = 2(kv2 = vi) bo" = yoo wiwo’ 
+y¥i(1 — Wo) Wo’ +yo(1 — Wo )pr’ 
and yields Y2(x) in a similar manner. Thus, 
ae é& 
u(x) =2 cos x-+—— sin* x +———- 
1+e (1+ .)- 
X[E— 1 cos x+ 3i6 cos 3x 


— 59, cos 5xv+sin *v]+---, (63) 
1 € is ¢€ i3 6 é 
ran ares rare armas meee 
I+e (1+e)? 16(1+6)? 16 (1+¢€)! (64) 


It is interesting to note that assuming e1 and 
expanding in (63) and (64) each coefficient in 
powers of e, we obtain the same series for u and 
vy as in Part B. 

Table I contains, in two columns, values of v 


TABLE I, 
v (according to 
€ vy (from series (64)) Van der Pol) 
0.1 0.998 ? 
0.33 98 0.99 
1 .93 0.90 
2 17 0.78 
5.6 A3 0.38(?) 
8 35 0.39 
10 .30 0.31 
11.4 ha 0.24 
99 .037 - 


for various values of e, calculated from the first 
four terms in (64), also from the experimental 
curves obtained by Van der Pol (for the latter 
1 am indebted to my Moore School colleague 
Professor C. N. Weygandt, Jr.). The agreement 
between the two columns is quite good. The 
first four terms in our series (64) give further 


ve~3.75 (ex), 


which again agrées well with Liénard’s theoretical 
results: ve~3.89---. Moreover, the same terms 
in (64) gives dvy/de<0, which makes plausible 
the following answer to a question raised by 
Liénard. As the parameter ¢ in Van der Pol’s 
equation increases from 0 to ~, the frequency » 
of its periodic solution decreases from 1 to 0, 
hence, the period T increases from 27 to ~. 

The foregoing method is offered for the purpose 
of computations. | wish to thank my pupil, Mr. 
J. D. Chapline, Jr., for checking the above 
numerical results for u(x) and ». 
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“Photographic Plates for Use in 
Spectroscopy and Astronomy” 


FIFTH EDITION 


new, fully revised edition of this well-known publication is 
A available. In it are listed the characteristics of the special plates for 
scientific work which are made in the Kodak Research Laboratories. There 
are over one hundred kinds of these plates, distributed among seven basic 
types of emulsion, each carefully prepared to provide desirable combinations 
of speed, contrast, granularity, and resolving power. They are sensitized in 
various ways to permit photography in spectral regions ranging from the 
short-wave-length ultraviolet to the infrared at 12,000A. 


A copy of the new edition of the booklet will be sent free upon request. 


EASTMAN KODAK COMPANY 


Research Laboratories Rochester 4, N. Y. 
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Arnold Engineering Company. 
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uve to The Arnold Engineering 
Company to continue with the same 


high devotion, energy and skill to 


Gehan mmelel products for the war effort. 
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Recent Applications of Physics 








Prepared by CLARK GOODMAN, Associate Editor 
Material to be included in this section shou!d be submitted to Dr. 


Clark Goodman, Massachusetts Institute of Technology, Cambridge 39, 
Massachusetts 


Electrostatic Dust 
Precipitator 


The cleaning action upon 
gases of a corona discharge 
has been known for more 
than a century. The Cottrell precipitator has long been 
recognized as the outstanding method of removing the 
very fine particles from gases, especially at elevated tem- 
peratures. The more general use of this technique for the 
removal of undesirable mist, dust, and pollen has been 
delayed for a number of reasons. The cost and maintenance 
are high as compared to other types of cleaning equipment. 
Considerable space and appreciable currents at d.c. voltages 
of from 30,000 to 100,000 are required. In addition, this 
type of precipitator generates so much ozone that the clean 
air, although free of dust, is too irritating to the nose and 
throat to be used for ventilation. Recently Westinghouse 
engineers have been able to overcome most of these objec- 
tions in an electrostatic unit known as the Precipitron. By 
careful design a compact unit small enough to be placed in 
an ordinary window has been developed. The generation Of 
ozone by a given discharge current depends on the geometry 
of the electrodes. By selecting a favorable design of elec- 
trode the generation of ozone per unit of current has been 
reduced to the point where its concentration in the air is 
comparable with that found in outdoor air in bright sunlight. 

The Precipitron has been utilized to remove oil mist 
from the atmosphere surrounding high speed grinders 
cooled by cutting oil. The expensive oil is recovered for 
re-use. Precipitrons have also been used to remove dust 
from the large volume of warm air required for dehydration 
of milk. The removal of dust from blast-furnace gases is 
sufficiently complete that these gases can be used as fuel 
for engines. Persons suffering from hay fever and certain 
forms of asthma are greatly relieved by the electrostatic 
removal of pollens from the air. Among the more promising 
postwar applications of electrostatic precipitation are the 
removal of road dust at the air intake of automobile engines 
and as a much needed supplement to present railway air- 
conditioning units. 
Mighty Midget Motor According to General Elec- 
tric engineers, the surface 
speed required to produce a high grade finish economically 
by grinding should be about the same whether grinding a 

large or a small hole. 
For 1-inch diameter 
or larger, medium-size 
wheels rotating at 
moderate speeds give 
the requisite surface 
speed. For small holes, 
many less than } inch 
in diameter, a_ tiny 
grinding wheel rotat- 
ing at a high speed is 


needed. A new G.E. motor has been developed for this pur- 
pose. It weighs but seven pounds, yet is rated 3 horsepower, 
With normal voltage applied at 2000 cycles, the midget 
motor reaches its full speed of 120,000 r.p.m. in less than a 
second. Tests have shown that these motors can deliver the 
rated power at this speed for more than eight hours of con- 
tinuous operation. The amount of material used in the 
motor is so small, with a consequent reduction in radiating 
surface, that water cooling is necessary, one-half a gallon per 
minute being used. For contrast, the small motor is com- 
pared in the accompanying photograph with the standard 
3-hp, 1800-r.p.m. motor which weighs 105 pounds. 


Modern three-dimensional 
photography has antiquated 
the old stereoscope that provided so much parlor entertain- 


Vectographic Views 


ment in the good old days. No longer is only one person 
at a time permitted to view a picture in three dimensions. 
Instead, an entire class or a large audience can see the 
same three-dimensional photograph or drawing when pro- 
jected on a screen and viewed through special polarizing 
spectacles. The techniques that make this possible were 
developed by the Polaroid Corporation. The products are 
called vectographs. They can be equally well used as prints 
or lantern slides. Both are filling urgent war needs in 
teaching military navigation and in aerial reconnaissance. 
In the training program a textbook containing vectographic 
prints is used in place of the conventional two-dimensional 
diagrams which attempt to represent spatial relationships. 
Celestial navigation vectograph slides, projected by stand- 
ard projectors, are so strikingly realistic that students feel 
they are looking at precise wire models of the heavens with 
relative positions of the stars and the earth immediately 
apparent. By revealing land details and the positions of 
enemy installations, vectographs enable key men in combat 
to know the lay of the land before they engage in offensive 
action against objective areas. Vectographs also offer the 
usual advantages of three-dimensional studies of bombing 
targets both before and after the attack. 

Vectographs are prepared from the two views or stereo- 
grams, required for all three-dimensional pictures, by a 
simple imbibition-transfer process. Prints of the stereo- 
scopic negatives—one normal and one reversed—are made 
on Eastman Wash-Off Relief Film. After development the 
picture images appear in relief on the gelatin surfaces. 
These reliefs, while still wet, are hinged together with their 
emulsion sides face to face and soaked briefly in a printing 
solution. A sheet of vectograph film, resembling a clear 
acetate sheet, is placed between them, and the sandwich 
is run through a wringer. The reliefs are then peeled from 
the vectograph film, which, as a result of these operations, 
has stereoscopic prints of the object on its two surfaces. 
Polarizing spectacles, in conjunction with the polarizing 
element in the vectograph print, permit the two eyes to 
see mutually exclusive images which the brain interprets 
as the normal seeing condition. As many as 2000 vecto- 
graphs can be made from a single set of relief films. The 
dried vectograph film may be mounted directly as a lantern 
slide or used as a print by painting the front surface with 
clear lacquer and the back surface with aluminum lacquer. 





